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PREFACE 


This report discusses work performed under NASA Contract 
NAS8-32234 for the George C. Marshall Space Flight Center. The 
cognizant technical officer was Mr. Frank Bugg, to whom we express 
our appreciation for his valuable assistance in obtaining the 
propellant for the dynamic property characterization testing. 
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ABSTRACT 


The results of a research study to characterize the dynamic 
response properties of the Space Shuttle SRM (TP-H1148) propellant 
and to establish the expected limits of propellant variability is 
presented in this final report. 

Dynamic shear modulus tests were conducted on six production 
batches of TP-H1148 at various static and dynamic strain levels over 
the temperature range from 40°F to 90°F. 

A heat conduction analysis and dynamic response analysis of 
the Shuttle SRM were also conducted. 

The dynamic test results show significant dependence on static 
and dynamic strain levels and considerable batch-to-batch and within- 
batch variability. However, the results of the SRM dynamic response 
analyses clearly demonstrate that the stiffness of the propellant has 
no consequential effect on the overall SRM dynamic response. Only the 
mass of the propellant needs to be considered in the dynamic analysis 
of the Shuttle SRM. 

Based on the results of the SRM dynamic analysis it is concluded 
that no special tests are required for quality control of the SRM 
propellant dynamic response properties. 
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DYNAMIC CHARACTERIZATION AND ANALYSIS OF 
SPACE SHUTTLE SRM SOLID PROPELLANT 


I. INTRODUCTION 


The program objectives were to perform a dynamic shear modulus 
characterization of TP-H1148 (Space Shuttle) propellant and to investi- 
gate the sensitivity of the Shuttle Solid Rocket Motor (SRM) dynamic 
response to propellant modulus variability associated with batch-to- 
batch and within-batch casting variations and temperature variations in 
the SRM. 

Specific tasks included: 

ASK 1: DYNAMIC SHEAR MODULUS CHARACTERIZATION 

A. Examine and analyze existing dynamic 
response test data for TP-H1148 
propellant and other similar solid 
propellants. 

B. Experimentally determine dynamic shear 
response properties of TP-H1148 pro- 
pellant as a function of frequency, 
temperature and strain level. 

C. Perform statistical analyses of dynamic 
shear data obtained to establish 
probability distributions associated 
with batch-to-batch and within-batch 
variability. 

D. Extend the SRM propellant dynamic 
response model to include the effects 
of propellant variability found in 

Task 1-B, and strain level, if necessary. 
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TASK 2: SRM DYNAMIC RESPONSE ANALYSES 

Examine the sensitivity of the SRM frequency 
and mode shape response to the propellant 
variability expected based on the results of 
Task 1 and to modulus gradients throughout the 
grain web associated with a nonuniform tempera- 
ture distribution and batch-to-batch propellant 
variability. 


TASK 3: IDENTIC . TEST REQUIREMENTS 

Based on the results of Tasks 1 and 2, recommend 
the tests required to characterize the dynamic 
response of each propellant casting and each SRM 
segment. 


The analysis effort and test data interpretation were carried out by 
W. L. Hufferd & Associates, and the characterization of TP-H1148 
propellant, supplied by Thiokol Corporation, Wasatch Division, was 
conducted by the Chemical Systems Division (CSD) of United Technologies 
Corporation under a subcontract from W. L. Hufferd & Associates. 

This final report presents the details of the rest *ch study. All 
technical objectives have been accomplished. The complete final report 
from CSD on the dynamic characterization of TP-H1148 propellant is 
included as Appendix A to this report. 
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II. SUMMARY OF ACCOMPLISHMENTS 

All technical requirements and objectives of the research study 
have been met. 

2.1 TASK 1: DYNAMIC SHEAR MODULUS CHARACTER I CAT I ON 

Propellant characterization tests were conducted by CSD on material 
from one 600 galon batch (TP-H1 1 48-9) and five other batches (TP-H1148- 
6,-7,-8,-9970096 and -9970115). The dynamic shear modulus of the six 
SRM solid propellant batches was measured at various static and dynamic 
strain levels. Constant elongation rate properties were also measured 
at 2 in/min at C°F, 40°F, 70°F, 90°F and 120°F, and the viscoelastic 
relaxation modulus was measured at 40°F, 70°F and 90°F. 

The test data show very small batch-to-batch propellant variability 
and excellent sample-to-sample reproducibility. Dynamic propellant 
modulus values were sensitive to the static compressive strain level as 
well as dynamic strain. The test strain levels covered 0.001% to 3% 
and a .orrespondingly large modulus range. This strain range should be 
applicable to the shuttle SRM dynamic loading conditions [1-5]. Specific 
modulus values for dynamic analyses can be selected from the data 
pre-cnted at the appropriate frequency and strain levels and used with 
i NASTRAN analysis. 

Additional statistical analyses of available test data for similar 
propellants and for TP-H1148 characterized by Thiokol /Wasatch Division 
r 5- 1 l] were performed to establish expected propellant variability. 
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Finally, the SRM propellant dynamic response model developed in 
[4] was extended to include the effects of propellant variability. 

2.2 TASK 2: SRM DYNAMIC RESPONSE ANALYSES 

Dynamic response analyses of a typical SRM segment were performed 
to examine the sensitivity of the SRM frequency and mode shape responses 
to propellant variability and to modulus variations associated with a 
non-uniform temperature gradient likely to exist in the SRM at the time 
of launch. It was found that the SRM dynamic response is insensitive to 
propellant variability and that no significant temperature gradients 
would exist through the SRM propellant web at the time of launch. The 
propellant provides most of the mass, but the case provides almost all 
of the stiffness due to the large differences between the propellant 
and case moduli at all temperatures, frequency and strain levels of 
interest in the Shuttle SRM. 

2.3 TASK 3: IDENTIFY TEST REQUIREMENTS 

Based on the insensitivity of the SRM dynamic response to propellant 
variability it is concluded that no special quality assurance testing 
is required for the dynamic response properties of TP-H1148 propellant. 

The normal quality control tests for ignition, burn rate and static 
mechanical properties are sufficient to indicate anomalous batch behavior. 
It must be emphasized, however, that anomalous behavior of one or even 
several castings of propellant, will still have an insignificant effect 
on the shuttle SRM dynamic response. 



III. DYNAMIC RESPONSE OF RELATED PROPELLANTS 


In addition to characterizing TP-H1 1 48 propellant during this 
program existing published data on TP-H1148 and similar related propel- 
lants were reviewed and analyzed to provide a broader data base for 
subsequent statistical analyses. Dynamic mechanical property tests 
have been conducted at the University of Utah [7] on the inert UTI-610 
PBAN propellant used in the NASA/Langley dynamic model tests [5]. 

Thiokol /Wasatch has characterized live and inert TP-H1123 PBAN 
propellant [6], which closely resembles the SRM propellant, as well as 
TP-H1148 propellant [10,11]. The Space Division of Rockwell Inter- 
national [8,9] has also conducted dynamic tests in support of one- 
quarter scale model tests; however, these data are difficult to analyze 
and draw conclusions from since non-standard tests were used and the 
test procedures are not well documented. 

3.1 UNIVERSITY OF UTAH T ESTS 

Characterization tests were conducted on the inert solid propellant 
cast in the 1/8 scale models of the space shuttle SRB dynamically tested 
at the NASA/Langley Research Center [5]. Static and dynamic tests were 
conducted to evaluate the interconvertibility of static and dynamic 
response properties according to linear viscoelasticity theory. These 
tests were carried out at several temperatures and the time (frequency)- 
temperature superposition principle used to extehd the range of 
applicability of the data. Some tests were also conducted to evaluate 
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the effects of moisture (i.e., relative humidity) on propellant response 
and also the effects of static strain level. 

The propellant, designated UTI-610, is an inert Polybutadiene 
Acryl o-Nitrile (PBAN) propellant with ammonium sulfate and sodium 
chloride replacing the normal oxidizer, ammonium perchlorate. The 
propellant samples and models were manufactured by the Chemical Systems 
Division of United Technologies Corporation, Sunnyvale, California. 

Propellant for the experimental characterization program was 
received from Nasa/Langeley Research Center in sealed one-half gallon 
milk cartons. Specimens 1.27 cm by 1.27 cm by 10.2 cm (0.5 in by 0.5 in 
by 4.0 in) were machined and wooden tabs bonded to the enc.s for the static 
testing, i.e., constant strain-rate and stress relaxation tests. Thin 
slabs, .1 cm by .3 cm by .2 cm (0.04 in by 0.12 in by 0.8 in) were used 
for the dynamic tension tests and slabs .15 cm by .15 cm by .3 cm 
(0.06 in by 0.06 in by 0.12 in) were used for the dynamic shear tests. 

Test specimen configurations are shown in Figure 1. 

A minimum of three and usually four replicate tests were carried 
out for each condition. 

The test methods employed conformed to standard industry-wide 
practices [12-14]. An Instron Universal testing machine (Figure 2) 
with Missimer temperature conditioning box was used for static testing 
and a Rheovibron dynamic tester (Figure 3) was used for the dynamic 
tests. 

Static and dynamic tests were conducted at -18°C (0°F), 4°C (40°F), 
25°C (77°F) and 49°C (120°F) and time- temperature superposition used 



WOODEN TABS (TYP) 



PROPELLANT 


STATIC TEST SPECIMEN 




DYNAMIC SHEAR TEST SPECIMEN 


FIGURE 1 - Test Specimen Configurations 
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FIGURE 3 - Schematic of Rheovlbron Dynamic Tester 
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to evaluate the response properties at other temperatures and 
frequencies.^ 

Constant strain rate tests were conducted at the strain rates 
and temperatures shown in Table I. True stress-strain response is 
shown in Figure 4. 

Stress relaxation tests were conducted according to Table II. 
Relaxation modulus as a function of log time and strain level is shown 
in Figure 5. With the exception of the one-percent strain level tests 
at 4°C the average values of the one and five-percent strain level test 
data were within the + 3a data band associated with the three-percent 
strain level tests. 

Figure 6 shows the time-temperature shift factor versus tempera- 
ture as determined from stress-relaxation tests, dyr .mic tension and 


It has been widely found that temperature has the effect of expand 
ing or contracting the time scale of response viscoelastic materials 
and that an equivalence between time or frequency and temperature exists 
Time and frequency are roughly the inverse of one another so that 
short time or high frequency response at one temperature corresponds to 
lower temperature. The converse holds true at higher temperatures. 

Thus, by obtaining propellant response at several temperatures, a 
time-temperature shift function, ay, relating the equivalence of tir.e 
or frequency and temperature can be experimentally determined by 
horizontally "shifting" the test data so that it superimposes to form, a 
single curve at some given reference temperature, usually 21 to 25°C. 

The resulting curve is known as the "master" response curve and is 
expressed in terms of temperature-reduced time, t/ay, in the case of the 
master-reduced frequency, way, in the case of the master dynamic moduli 
curves (See i.e., [12-15]). 

Thus, conventional analyses remain valid in terms of the master 
relaxation and dynamic moduli if time, t, and frequency, w, are replaced 
by temperature-reduced time, t/ay, and frequency-reduced time, way, 
respectively. 1 
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TABLE I 

CONSTANT STRAIN RATE TESTS 


N \ v Temperature 
StralrNs. (°C) 

Rate i 
(min)" 1 

-18 

4 

25 

49 

5 



X 

X 

0.5 

X 

X 

X 

X 

0.05 


X 

X 


0.005 

X 


X 










TRUE STRESS (kN/m 2 ) 


2 - 


3850 


3300 


2’50 


2200 


1650 


1100 


550 


0 
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TABLE II 

STRESS RELAXATION TESTS 








FIGURE 5 - Stress Relaxation Modulus Versus Log Time 


LOG 
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TEMPERATURE (°C) 


FIGURE 6- Time-Temperature Shift Factor 
Versus Temperature 
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dynamic shear tests. Three-sigma limits are also shown for the stress- 
relaxation determined shift factors. 

Thermorheologically simple material behavior requires time- 
temperature shift factors to be identical, irrespective of how they 
were determined. However, it may be noted in Figure 6 that the shift 
factors determined from dynamic shear tests do not agree with those 
determined from stress-relaxation or dynamic tension tests, which are 
in reasonable agreement with each other. The reason for this discrep- 
ancy is not known, although it is probably associated with some 
deficiency in the dynamic shear test procedure. The shift factors 
associated with the stress relaxation tests were used to superpose the 
relaxation modulus and dynamic tension test results whereas the shift 
factors obtained from dynamic shear tests were used to develop a master 
dynamic shear modulus curve. 

Figure 7 presents the master relaxation modulus as a function of 
temperature-reduced time, t/a^, a + 25°C. Figure 7 also shows the 
master relaxation modulus obtained from the constant strain rate test 
data shown in Figure 4 using the relation [15], 


W t/a T> * 


da(t/a T ) I 


de 


:=Rt/a, 


where R is the rate of straining. It may be observed that the modulus 
obtained from'the constant strain rate tests is substantially higher 
than that obtained from stress relaxation tests, particularly at short 
times or low temperatures. This nonlinear behavior is typically observed 
of solid propellants. 




FIGURE 7 - Master Stress Relaxation 
Modulus Versus Temperature Reduced Time 
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Dynamic tension tests were conducted at 3.5, 11, 35, and 110 Hz 
and 0.1, 0.5, 1.0, 2.5, and 5.0 percent static strain levels at 18°C 
(0°F) , 4°C (40°F), 25°C (77°F) and 49°C (120°F). Dynamic shear tests 
were conducted at the same frequencies and temperatures at pre-imposed 
static strain levels of 0.5, 1.0, 2.5 and 5.0 percent. 

A fixed half-amplitude dynamic displacement of 50 m was used for 
all dynamic tests. The length of the dynamic specimens for the 0.1 
percent strain level tests was 5 cm (2 in) rather than the 2 cm 
mentioned previously. For this test the dynamic strain level was equal 
to the imposed static strain of 0.1 percent. In all other dynamic 
tests the dynamic strain level varied between 0.22 and 0.25 percent. 

Dynamic tension results are tabulated in Tables Ill-a through Ill-d 
and dynamic shear results are tabulated in Tables IV-a through IV-d. 

The notation in Tables III and IV is derived from the following repre- 
sentations for complex modulus: 


E*(w) = E' + IE" 
G*(u>) =6’ + iG" 



| G | = V(S') 2 + (G‘) 2 


It may be observed from Tables III and IV that a gradual decline 
or softening of the dynamic stiffness occurs with increasing static 
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TABLE Ill-a 

DYNAMIC TENSION TEST RESULTS 
AT A FREQUENCY OF 3.5 H* 


STATIC 

STRAIN 

(X) 

TEMPERATURE 

(°C) 

tan 5 

|E| , 

(MH/m) 

E ' 2 
(MN/in ) 

t #i 

[m/m 2 

0.1 

-18 

0.32 

313 

299 

95.1 


4 

0.47 

65.3 

59.0 

27.8 


25 

0.45 

16.5 

15.0 

6.74 


49 

0.38 

3.92 

3.68 

1.36 

0.5 

-18 

0.32 

179 

170 

54.4 


4 

0.45 

45.6 

41.5 

18.8 


25 

0.43 

13.5 

12.4 

5.3 


49 

0.32 

5.91 

5.58 

1.82 

1.0 

-18 

0.33 

175 

166 

53.7 


4 

0.45 

43.8 

40 

17.9 


25 

0.43 

12.6 

11.6 

4.6/ 


49 

0.30 

5.42 

5.17 

1.58 

2.5 

-18 

0.34 

166 

158 

52.3 


4 

0.45 

39.9 

36.6 

16.3 


25 

0.42 

11.6 

10.53 

4.34 


49 

0.30 

4.68 

4.48 

1.35 

5.0 

-18 

0.35 

138 

131 

42.3 


4 

0.43 

34.8 

32.1 

13.4 


25 

0.39 

10.9 

9.42 

3.11 


49 

0.26 

4.0 

3.85 

1.09 
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TABLE Ill-b 

DYNAMIC TENSION TEST RESULTS 
AT A FREQUENCY OF 11 Hz 


STATIC 

STRAIN 

(%) 

TEMPERATURE 

(°C) 

tan 5 

Imi/m 2 ) 

E ‘ 2 
(MN/nr) 

E" 

(MN/m 2 ) 

0.1 

-18 

0.34 

333 

315 

107 


4 

0.50 

70.8 

63.3 

31.5 


25 

0.47 

19.7 

17.8 

8.35 


49 

0.42 

6.75 

6.17 

2.72 

0.5 

-18 

0.33 

207 

197 

63.3 


4 

0.48 

53.8 

48.5 

23.5 


25 

0.48 

15.0 

13.6 

6.49 


49 

0.38 

5.99 

5.6 

2.11 

1.0 

-18 

0.33 

214 

204 

66.4 


4 

0.49 

52.5 

47.2 

23.2 


25 

0.47 

14.0 

12.6 

6.0 


49 

0.37 

6.22 

5.85 

2.13 

2.5 

-18 

0.32 

213 

204 

62.6 


4 

0.48 

55.8 

50.33 

24.2 


25 

0.47 

14.7 

13.2 

6.32 


49 

0.36 

6.48 

6.08 

2.24 

5.0 

-18 

0.33 

185 

177 

52.5 


4 

0.46 

48.5 

44.0 

20.5 


25 

0.44 

13.8 

12.6 

5.64 


49 

0.35 

5.88 

5.54 

1.97 
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TABLE III-c 

DYNAMIC TENSION TEST RESULTS 
AT A FREQUENCY OF 35 Hz 


STATIC 

STRAIN 

(*) 

TEMPERATURE 

(°C) 

tan 5 

Iwi/m 2 ) 

E ' 2 
(MN/nr) 

E M 

(MN/m 2 ) 

0.1 

-18 

0.42 

405 

374 

153 


4 

0.68 

71.2 

58.9 

39.7 


25 

1.08 

25.0 

17.2 

18.0 


49 

0.64 

8.8 

7.43 

4.72 

0.5 

-18 

0.26 

269 

260 

67.3 


4 

0.55 

75.2 

65.9 

36.2 


25 

0.62 

21.8 

18.5 

11.4 


49 

0.47 

8.37 

7.56 

3.57 

1.0 

-18 

0.26 

276 

267 

69.3 


4 

0.55 

75.2 

65.9 

36.2 


25 

0.62 

21.8 

18.5 

11.4 


49 

0.47 

8.37 

7.56 

3.57 

2.5 

-18 

0.26 

263 

254 

66.1 


4 

0.52 

76.3 

68 

34.7 


25 

0.58 

21.6 

18.8 

10.7 


49 

0.46 

8.66 

7.85 

3.64 

5.0 

-18 

0.27 

243 

234 

61.8 


4 

0.48 

71 

64 

30.8 


25 

0.50 

18.6 

16.7 

8.24 


49 

0.45 

7.92 

7.23 

3.24 
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TABLE Ill-d 

DYNAMIC TENSION TEST RESULTS 
AT A FREQUENCY OF 110 Hz 


STATIC 

STRAIN 

(2) 

TEMPERATURE 

(°C) 

tan & 

|EI , 
(MN/nr) 

V 2 

(MN/nr) 

E M 

(MN/m 2 ) 

0.1 

-18 

0.375 

512 

482 

171 


4 

1.1 

120 

87.2 

78.1 


25 

1.24 

33.5 

27.1 

24.2 


49 

1.25 

12.3 

8.06 

8.94 

0.5 

-18 

0.27 

316 

304 

81.4 


4 

0.53 

102 

90.4 

48.0 


25 

0.59 

30.2 

26.0 

15.4 


49 

0.54 

12.2 

10.8 

5.81 

1.0 

-18 

0.26 

328 

308 

80.7 


4 

0.53 

105 

92.6 

49.2 


r- 

0.58 

31.1 

26.8 

15.8 


49 

0.53 

12.2 

10.7 

5.72 

2.5 

-18 

0.25 

324 

314 

75.8 


4 

0.50 

108 

96.9 

47.8 


25 

0.56 

31.6 

27.6 

15.5 


49 

0.51 

12.3 

11.3 

5.59 

5.0 

-18 

0.26 

299 

289 

73.4 

4 

0.48 

103 

92.0 

44.4 


25 

0.52 

30.3 

2*. a 

14.0 


49 

0.48 

11.7 

10.6 

5.12 
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TABLE IV-a 

DYNAMIC SHEAR TEST RESULTS 
AT A FREQUENCY OF 3.5 H z 


STATIC 

STRAIN 

(%) 

TEMPERATURE 

(°C) 

tan <5 

ImN/id 2 ) 

G ’ ? 
(MN/nT) 

G' 1 

(MN/m 2 ) 

0.5 

-18 

0.052 

11.2 

11.2 

0.59 


4 

0.174 

8.06 

7.94 

1.33 


25 

0.29 

4.33 

4.16 

1.18 


49 

0.31 

1.89 

1.81 

0.55 

1.0 

-18 

0.061 

11.0 

11.0 

0.67 


4 

0.21 

7.35 

7.16 

1.49 


25 

0.33 

3.44 

3.27 

1.05 


49 

0.31 

1.71 

1.64 

0.51 

2.5 

-18 

0.061 

11.1 

11.1 

0.68 


4 

0.21 

7.29 

7.13 

1.48 


25 

0.33 

3.26 

3.12 

1.01 


49 

0.32 

1.59 

1.52 

0.48 

5.0 

-18 

0.059 

11.1 

11.1 

0.66 


4 

0.205 

7.35 

7.2 

1.45 


25 

0.33 

3.26 

3.1 

1.01 


49 

0.32 

1.54 

1.47 

0.47 
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TABLE IV-b 

DYNAMIC SHEAR TEST RESULTS AT 
A FREQUENCY OF 11 Hz 


STATIC 

STRAIN 

.(*) 

TEMPERATURE 

CO 

tan 6 

ImiI/hi 2 ) 

G 1 

(MN/m 2 ) 

G" 

(MN/m 2 ) 

0.5 

-18 

0.023 

12.63 

12.6 

0.30 


4 

0.10 

10.5 

10.4 

1.08 


25 

0.24 

6.5 

6.32 

1.51 


49 

0.31 

3.42 

3.28 

0.99 

1.0 

- 1 8 

0.023 

12.8 

12.8 

0.31 


4 

0.1 

10.4 

10.3 

1.14 


25 

0.27 

6.0 

5.78 

1.55 


49 

0.33 

2.99 

2.67 

0.95 

2.5 

-18 

0.023 

12.8 

12.8 

0.31 


4 

0.11 

10.5 

10.4 

1.11 


25 

0.27 

5.86 

5.64 

1.59 


49 

0.34 

2.82 

2.67 

0.91 

5.0 

-18 

0.023 

12.8 

12.8 

0.30 


4 

0.093 

10.5 

10.4 

1.02 


25 

0.28 

5.55 

5.33 

1.53 


49 

0.34 

2.58 

2.43 

0.85 
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TABLE IV-c 

DYNAMIC SHEAR TEST RESULTS 
AT A FREQUENCY OF 35 Hz 



TEMPERATURE 

(°c) 

tan 6 

0.5 

-18 

0.036 


4 

0.15 


25 

0.31 


49 

0.41 

1.0 

-18 

0.045 


4 

0.22 


25 

0.41 


49 

0.42 

2.5 

-18 

0.05 


4 

0.23 


25 

0.44 


49 

1 . 42 

5.0 

-18 

0.072 


4 

0.28 


25 

0.45 


49 

0.42 


ImN/ih 2 ) 

G 1 

(MN/m 2 ) 

G" 2 
(MN/nr) 

14 67 

14.67 

0.58 

11.5 

11.4 

1.82 

6.26 

5.97 

1.85 

2.46 

2.29 

0.91 

14.3 

14.3 

0.68 

9.83 

9.58 

2.18 

3.8 

3.52 

1.43 

1.59 

1.47 

0.61 

14.0 

13.9 

0.75 

9.13 

8.85 

2.16 

3.28 

3.01 

1.28 

1.2° 

1.13 

0.47 

13.2 

13.2 

1.03 

7.8 

7.5 

2.09 

2.97 

2.72 

1.23 

0.99 

0.91 

0.37 
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TABLE IV-d 

DYNAMIC SHEAR TEST RESULTS AT A 
FREQUENCY OF 110 Hz 


STATIC 

STRAIN 

(*) 

TEMPERATURE 

ro 

tan 6 

ImN/bi 2 ) 

G ' 2 
(MN/nr) 

G* 1 

(MN/m 2 ) 

0.5 

-18 

0.036 

22.0 

22.0 

0.34 


4 

0.12 

18.0 

17.9 

2.13 


25 

0.29 

10.9 

10.4 

2.98 


49 

0.41 

4.99 

4.65 

1.91 

1.0 

-18 

0.046 

22.2 

22.2 

1.08 


4 

0.14 

17.6 

17.4 

2.50 


25 

0.34 

9.77 

9.25 

3.08 


49 

0.43 

4.48 

4.12 

1.77 

2.5 

•18 

0.046 

22 4 

22.4 

1.12 


4 

0 15 

17.2 

U.O 

2.58 


25 

0.35 

9.33 

8.8 

3.09 


49 

0.45 

4.01 

3.65 

1.66 

5.0 

-18 

0.049 

22.0 

22.0 

1.14 


4 

0.16 

16.8 

16.6 

2.66 


25 

0.37 

8.42 

7.96 

2.94 


49 

0.46 

3.07 

2.79 

1.30 



-27- 


strain level. The maximum dynamic tensile modulus occurs for a static 
strain of 0.1 percent where the dynamic strain is also equal to this 
value. Typical of live propellant behavior, dynamic moduli increase 
lor decreasing static strain and/or dynamic strain levels (see, e.g., 
[4,16,17]). 

Figure 8 presents the master dynamic tensile and shear moduli 
versus temperature reduced frequency at a reference temperature of 
25°C. The dynamic tensile modulus, calculated from stress relaxation 
test data and constant strain rate test data is also shown for comparison. 
It may be noted that the dynamic moduli calculated in this fashion do 
not agree with the measured dynamic data. This is typical of solid 
propellant behavior. The master dynamic shear modulus curve shewn in 
Figure 9 was obtained using the frequency-temperature shift factors 
developed from the dynamic shear tests and previously presented in 
Figure 6. Although considerable data scatter exists, it ray be seen 
that the master dynamic shear modulus is approximately one-third of 
the value o f the measured tensile modulus, indicating nearly incompres- 
sible material behavior. From bulk compressibility tests conducted by 

Thiokol/Wasatch, the hulk modulus was experimentally determined to be 

2 

267,000 psi (1.84 MN/m ) implying a Poisson's ratio in the range of 
0.4991 to 0.4999 over the frequency range of interest. 

The imaginary part of the dynamic tensile and shear modulus, E" and 
G" , are presented in Figure 9 as a function of temperature-reduced 
frequency. The variability and data scatter in G" indicate that the 
shear test results are not likely to be truly representative of the 
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propellant's dynamic stress behavior. The anomalous shear behavior is 
attributed to coupling between the shim-stock support and glue and the 
propellant (see Figure 1). 

3.2 THIOKOL/WASATCH TESTS 

Thiokol has conducted dynamic shear tests on live and inert 
TP-1123 propellant [6], which, as noted previously, is very similar to 
TP-H1148, and static ar.d dynamic characterization of TP-H1148 
propellant [10,11]. 

Tfiokol conducts dynamic shear tests using a modified Gottenberg 
disk test technique [12]. The test setup consistr of a disk of propel- 
lant, bonded to steel outer and center rings (Figure 10), and mechanically 
loaded at its center with an electromechanical vibrator at a constant 
displacement. Excitation force, disk center displacement and relative 
phase angle are measured from which dynamic shear moduli are calculated. 

Tests on live and inert TP-H1123 propellant were conducted at 
frequencies ranging from 10 to 50 Hz at 40°F, 70°F and 90°F. Analytical 
studies of the Shuttle SRM center segment and the Gottenberg disk test 
specimen were performed to determine appropriate static and dynamic 
displacements for the test. A Gottenberg disk static center-body 
deflection of 0.015 inch was determined to simulate the propellant 
strains resulting from the Worst case SRM thermal shrinkage, and a 
vibratory center-body displacement of 0.0006 inch peak-to-peak was 
required to simulate maximum propellant strains from a 0.5 g peak-to-peak, 
2 to 50 Hz, longitudinal, sinusoidal vibration input to the case. Since 
the above dynamic displacement was below the minimum amplitude limit of 
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FIGURE 10. Gottenberg Disk Test Specimen 
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the test apparatus, tests were conducted at dynamic displacement 
amplitudes ranging from 0.001 to 0.004 inch, and imposed center-body, 
static displacements ranging from 0.009 to 0.022 inches. 

Experimental test results as a function of dynamic displacement 
are presented in Table V for live TP-H1123 propellant >nd in Table VII 
for inert TP-H113 propellant (designated H-l 3) . Similar results are 
presented in Tables VII and VIII as a function of static displacement. 
Figure 11 presents the dynamic shear moduli, G' and G", versus tempera- 
ture reduced frequency for TP-H1123 and H-l 3 propellants. 

While the test data presented in Tables V through VIII indicate a 
significant influence of temperature and frequency, as expected; a 
Duncan Multiple Range Test indicated no statistical influence of static 
or dynamic displacement amplitude over the range studied. 

In addition, Thiokol also conducted Gottenberg disk tests on TP-H1123 
live propellant in a pressure vessel at ambient, 500 and 1000 psig 
confining pressures. Their results are presented in Table IX. These 
data are well within the statistical limits established by the ambient 
pressure tests previously presented. 

Thiokol has also conducted static [10] and dynamic [11] charac- 
terization tests of TP-H1148 space shuttle propellant. Dynamic shear 
tests were conducted on propellant from six mixes (9910-053, -057, 

-059, -060, -061, -062) at 40, 70 and 9C°F and frequencies ranging from 
5 to 50 Hz. The tests were conducted using a Gottenberg disk with a 
dynamic displacement of 0.002 inches. The test data are summarized in 
Tables X through XV, and analyzed subsequently. 
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TABLE V 

DYNAMIC SHEAR MODULUS OF LIVE TP-H 1123 PROPELLANT 
TEST DATA FOR VARIOUS DYNAMIC DISPLACEMENTS 


Dynamic 

Displacement 

(Inches) 

Frequency 

(Hz) 

Temperature 

(°F) 

tar.6 

, |G !, 

(psi) 

G' 

(psi) 

G" 

(psi) 

0.001 

50 

40 

.44 

3597 

3064 

1363 


40 


.52 

3014 

2670 

1398 


30 


.47 

2573 

2328 

1095 


20 


.38 

2168 

2029 

765 


15 


.32 

1949 

1856 

596 


10 


.25 

1757 

1705 

423 

0.002 

50 

40 

.53 

3294 

2909 

1545 



70 

.58 

1472 

1206 

704 



90 

.62 

801 

680 

423 


40 

40 

.52 

2856 

2532 

1322 



70 

.58 

1178 

1019 

590 



90 

.61 

660 

562 

345 


30 

40 

.47 

2426 

2196 

1030 



70 

.53 

974 

859 

457 



90 

.57 

331 

402 

262 


20 

40 

.38 

2017 

1888 

710 



70 

.41 

782 

722 

298 



90 

.47 

414 

374 

177 


15 

40 

.32 

1824 

1737 

558 



70 

.33 

680 

656 

220 



90 

.39 

366 

334 

130 


10 

40 

.25 

1632 

1584 

393 



70 

.24 

606 

590 

140 



90 

.28 

298 

287 

79.7 

0.0027 

50 

40 

.53 

3292 

2908 

1544 


40 


.52 

2852 

2592 

1319 


30 


.47 

2438 

2200 

1032 


20 


.38 

2026 

1896 

714 


15 


.32 

1828 

1740 

558 


10 


.24 

1655 

1608 

392 
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TABLE V (continued) 

DYNAMIC SHEAR MODULUS OF LIVE TP-H1123 PROPELLANT 
TEST DATA FOR VARIOUS DYNAMIC .DISPLACEMENTS 


Displacement Fr ft“?" c » Tan ?!" ture 

/ Inelme^ ' MZ > ' N 


0.003 50 70 

90 

40 70 

90 

30 70 

90 

20 70 

90 

15 70 

90 

10 70 

90 

0.004 50 70 

90 

40 70 

90 

30 70 

90 

20 70 

90 

15 70 

90 
70 
90 


tan6 

J G L 

(psi) 

G' 

(psi ) 

G" 

(psi) 

.58 

1289 

1116 

645 

.62 

793 

674 

417 

.58 

1089 

942 

544 

.61 

655 

559 

342 

.53* 

904 

797 

425 

.57 

526 

458 

259 

.42 

733 

676 

283 

.48 

397 

357 

173 

.34 

699 

618 

209 

.43 

322 

295 

128 

.24 

547 

532 

127 

.30 

277 

265 

80.3 

.57 

1315 

1140 

656 

.52 

777 

663 

406 

.56 

1106 

960 

542 

.61 

645 

552 

335 

.52 

908 

807 

416 

.47 

398 

360 

169 

C\1 

726 

669 

282 

.47 

398 

360 

169 

.35 

636 

601 

208 

.39 

343 

320 

125 

.27 

555 

539 

148 

.27 

288 

278 

79.0 


10 
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TABLE VI 

DYNAMIC SHEAR MODULI OF H-13 INERT PROPELLANT 
TEST DATA FOR VARIOUS DYNAMIC DISPLACEMENTS 


Dynamic 

Displacement 

(Inches) 

O.OOll 


0.002 


Frequency 

(Hz) 

Temperature 

(°F) 

tan6 

J 6 L 

(psi) 

G' 

(psi) 

G*‘ 

(psi) 

50 

40 

.52 

2504 

2220 

1153 


70 

.64 

998 

841 

538 


90 

.77 

710 

561 

430 

40 

40 

.51 

2197 

1953 

998 


70 

.63 

836 

708 

445 


90 

.77 

599 

473 

364 

30 

40 

.50 

1909 

1707 

849 


70 

.58 

673 

588 

340 


90 

.73 

494 

397 

290 

20 

40 

.46 

1642 

1486 

689 


70 

.47 

533 

482 

227 


90 

.63 

398 

335 

211 

15 

40 

.44 

1510 

1331 

604 


70 

.41 

46? 

428 

175 


90 

.55 

349 

304 

148 

10 

40 

.41 

1382 

1277 

519 


70 

.31 • 

400 

382 

118 


90 

.43 

305 

280 

119 

50 

40 

.52 

2292 

2030 

1060 


70 

.64 

995 

838 

536 


90 

.58 

616 

532 

310 

40 

40 

.51 

1981 

1764 

897 


70 

.63 

837 

708 

445 


90 

.56 

507 

443 

247 

30 

40 

.47 

1676 

1514 

714 


70 

.57 

680 

592 

336 


90 

.5 

407 

364 

182 
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TABLE VI (continued) 

DYNAMIC SHEAR MODULI OF H-13 INERT PROPELLANT 
TEST DATA FOR VARIOUS DYNAMIC DISPLACEMENTS 


Dynamic 

Displacement 

(Inches 

Frequency 

(Hz) 

Temperature 

(°F) 

tan 

G 

(psi) 

G' 

(psi ) 

G" 

(psi) 

0.002 

20 

40 

.43 

1391 

1277 

545 



70 

.47 

536 

484 

228 



90 

.44 

315 

288 

126 


15 

40 

.36 

1252 

1177 

422 



70 

.31 

396 

378 

116 



90 

.33 

232 

220 

72.5 

0.003 

50 

70 

.77 

946 

751 

575 


40 


.75 

796 

636 

478 


30 


.73 

648 

524 

382 


20 


.61 

506 

433 

262 


15 


.50 

439 

392 

197 


10 


' .37 

371 

349 

130 
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TABLE VII 

DYNAMIC SHEAR MODULI OF TP-H1123 LIVE PROPELLANT 
AS AFFECTED BY VARIOUS STATIC DISPLACEMENTS 


Static 


Displacement (inches) 

0.0098 

0.0177 

0.0217 

0.0197 

Frequency 





(Hz) 



G'(psi) 


50 

1153 

1178 

1163 

1153 

40 

982 

1003 

992 

982 

30 

828 

850 

839 

828 

20 

694 

723 

718 

694 

15 

627 

650 

650 

627 

10 

558 

588 

582 

558 




G"(psi) 


50 

656 

674 

663 

656 

40 

551 

565 

558 

551 

30 

421 

433 

427 

42: 

15 

197 

204 

204 

197 

10 

122 

128 

127 

122 




M(psl) 


50 

1327 

1357 

1339 

1327 

40 

1126 

1151 

1138 

1126 

30 

929 

954 

941 

923 

20 

745 

777 

771 

745 

15 

657 

681 

681 

657 

10 

571 

602 

596 

871 


Note: Tested at 70°F, ambient pressure and a dynamic displacement 

of 0.002 inch peak to peak. 
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TABLE VIII 

DYNAMIC SHEAR MODULI OF H-13 INERT PROPELLANT 
AS AFFECTED BY VARIOUS STATIC DISPLACEMENTS 


Static 

Displacement (inches) 0.0217 

Frequency 

.(Hz) 

50 1046 

40 89S 

30 765 

20 652 

15 588 

10 532 


50 602 

40 506 

30 393 

20 261 

15 195 

10 121 


50 1207 

40 1032 

30 860 

20 702 

15 619 

10 546 


0.0098 0.0086 0.0017 


G'Cpsi) 


986 

1076 

1001 

837 

925 

853 

706 

798 

722 

589 

675 

594 

523 

617 

535 

477 

554 

471 

(T 

(£§11 


557 

624 

568 

464 

524 

475 

359 

412 

368 

234 

271 

243 

173 

205 

177 

108 

128 

107 

G*(psi ) 


1132 

1244 

1151 

957 

1063 

976 

792 

898 

810 

634 

727 

642 

551 

650 

564 

489 

569 

483 


Note: Tested at 70°F and ambient pressure 



LOG G'M & G (cj) (psi) 



TABLE IX 

DYNAMIC SHEAR MODULI OF TP-H1123 LIVE 
PROPELLANT AS AFFECTED BY PRFSSURE 
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Note: Tested at 70°F and at a dynamic displacement of 0.0045 inch peak to peak. 
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TABLE X 

DYNAMIC SHEAR MODULUS, G', VERSUS FREQUENCY AT 90°F 


FREQUENCY (Hz) 


MIX 

(9910-) 5 

10 

15 

20 

30 

40 

50 

053 407 

525 

635 

742 

965 

1190 

1450 

455 

550 

670 

782 

1020 

1360 

1590 

328 

425 

521 

625 

830 

1030 

1270 

MEAN 297 

500 

609 

716 

938 

1160 

1440 

STD. DEV. 64 

66 

78 

82 

98 

118 

160 

057 298 

435 

545 

660 

880 

1080 

1310 

356 

440 

535 

635 

830 

1030 

1280 

278 

365 

460 

560 

749 

950 

1190 

MEAN 311 

413 

513 

618 

820 

1020 

1260 

STD. DEV. 41 

42 

46 

52 

66 

66 

62 

059 291 

390 

485 

587 

785 

1260 

1210 

291 

400 

480 

566 

754 

1100 

1180 

327 

420 

515 

630 

840 

948 

1290 

MEAN 303 

403 

493 

594 

793 

1100 

1230 

STD. DEV. 21 

15 

19 

33 

44 

156 

57 

060 416 

549 

654 

785 

1030 

1090 

1530 

340 

450 

567 

669 

890 

1110 

1340 

268 

355 

450 

550 

750 

1230 

1150 

MEAN 341 

451 

557 

668 

890 

1140 

1340 

STD. DEV. 74 

97 

102 

118 

140 

76 

190 

061 238 

400 

520 

628 

862 

1090 

1360 

360 

449 

555 

660 

883 

1110 

1400 

429 

532 

642 

760 

990 

1230 

1530 

MEAN 359 

460 

572 

683 

SI 2 

1140 

1430 

STD. DEV. 71 

67 

63 

69 

69 

76 

89 

062 357 

430 

535 

647 

860 

1080 

1360 

334 

437 

538 

648 

845 

1050 

1300 

325 

410 

500 

600 

800 

1000 

1220 

MEAN 339 

426 

524 

632 

835 

1040 

1290 

STD. DEV. 17 

14 

21 

27 

31 

40 

70 

GRAND MEAN 542 

442 

545 

652 

865 

1080 

1330 

STD. DEV. 55 

59 

66 

73 

88 

101 

128 
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TABLE XI 

DYNAMIC SHEAR MODULUS, G\ VERSUS FREQUENCY AT 90°F 


FREQUENCY (Hz) 
MIX 


(9910-) 

5 

10 

15 

20 

30 

40 

053 

243 

321 

400 

500 

690 

770 


240 

320 

417 

515 

718 

925 


177 

210 

294 

375 

555 

730 

MEAN 

220 

284 

370 

463 

654 

842 

STD. DEV. 

37 

64 

67 

77 

87 

101 

057 

167 

235 

315 

409 

590 

785 


235 

280 

350 

430 

595 

760 


156 

195 

260 

340 

500 

670 

MEAN 

186 

237 

308 

393 

562 

738 

STD. DEV. 

43 

43 

45 

47 

53 

60 

059 

176 

219 

280 

360 

530 

700 


160 

217 

281 

343 

505 

692 


175 

225 

305 

395 

525 

760 

:;ean 

170 

220 

289 

366 

537 

717 

STD. DEV. 

9 

4 

14 

27 

35 

37 

060 

180 

275 

370 

465 

661 

860 


175 

230 

312 

405 

600 

797 


151 

200 

265 

350 

510 

682 

MEAN 

169 

235 

316 

407 

590 

780 

STD. DEV. 

16 

38 

53 

53 

76 

90 

061 

136 

215 

305 

400 

590 

770 


149 

205 

290 

385 

580 

777 


186 

262 

355 

455 

646 

850 

MEAN 

157 

227 

317 

413 

605 

799 

STD. DEV. 

26 

30 

34 

37 

36 

44 

062 

165 

249 

317 

395 

563 

747 


165 

220 

300 

380 

552 

726 


156 

220 

290 

370 

540 

700 

MEAN 

162 

230 

302 

382 

552 

724 

STD. DEV. 

5 

17 

14 

13 

12 

24 

GRAND MEAN 

177 

239 

317 

404 

583 

967 

STD. DEV. 

31 

38 

44 

51 

62 

72 


50 


1100 

1170 

931 

1070 

123 

1000 

941 

858 

933 

71 

878 

878 

954 

903 

44 

1060 

1000 

851 

970 
108 

990 

1020 

1060 

1020 

35 

975 

920 

890 

928 

43 

971 
88 
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TABLE XII 

DYNAMIC SHEAR MODULUS, G' , VERSUS FREQUENCY AT 70°F 


FREQUENCY (Hz) 
MIX 


(9910 

-) 

5 

10 

15 

20 

30 

40 

50 

053 


597 

775 

910 

1060 

1350 

1640 

2000 



575 

7*3 

912 

1080 

1360 

1750 

1950 



680 

875 

1050 

1220 

1530 

1850 

2240 


MEAN 

617 

792 

960 

1120 

1410 

1710 

2060 


STD. DEV. 

55 

76 

78 

87 

101 

118 

155 

057 


570 

714 

850 

1000 

1270 

1570 

1920 



550 

750 

900 

1050 

1350 

1670 

2000 



634 

740 

890 

1050 

1350 

1670 

2020 


MEAN 

585 

735 

880 

1030 

1320 

1640 

1980 


STD. DEV. 

44 

19 

26 

29 

46 

58 

53 

059 


484 

650 

812 

950 

1250 

1550 

1940 



597 

740 

875 

1025 

1300 

1580 

1890 



650 

827 

1000 

1160 

1500 

1830 

2230 


MEAN 

577 

739 

896 

1040 

1350 

1650 

2020 


STD. DEV. 

85 

89 

96 

106 

132 

154 

184 

060 


616 

800 

950 

1100 

1420 

1740 

2130 



604 

790 

960 

1120 

1440 

1750 

2100 



544 

680 

830 

970 

1280 

1580 

1940 


MEAN 

588 

757 

913 

1068 

1380 

1690 

2060 


STD. DEV. 

39 

67 

72 

81 

87 

95 

102 

061 


564 

739 

900 

1080 

1410 

1750 

2200 



582 

700 

840 

990 

1280 

1580 

1950 



650 

900 

1080 

1260 

1600 

1940 

2350 


MEAN 

612 

780 

940 

1110 

1430 

1760 

2170 


STD. DEV. 

68 

106 

125 

137 

161 

180 

202 

062 


611 

750 

875 

1020 

1320 

1620 

1900 



546 

700 

860 

1020 

1320 

1640 

2020 



551 

705 

867 

1010 

1300 

1620 

2000 


MEAN 

569 

718 

867 

1020 

1310 

1630 

1970 


STD. DEV. 

36 

23 

8 

6 

12 

12 

64 


GRAND MEAN 

591 

753 

909 

1060 

1370 

1680 

2040 


STD. DEV. 

51 

65 

74 

83 

97 

109 

135 
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TABLE XIII 

DYNAMIC SHEAR MODULUS, G", VERSUS FREQUENCY AT 70°F 

FREQUENCY (Hz) 


MIX 


(9910-) 

5 

10 

15 

20 

30 

40 

50 

053 

350 

525 

660 

792 

1020 

1250 

1450 


347 

430 

560 

675 

925 

1180 

1460 


403 

514 

650 

775 

1040 

1280 

1580 

MEAN 

367 

490 

623 

747 

995 

1240 

1500 

STD. DEV. 

32 

52 

55 

63 

61 

51 

72 

057 

264 

367 

475 

600 

850 

1114 

1385 


225 

356 

475 

600 

850 

1125 

1400 


310 

400 

525 

640 

900 

1190 

1550 

MEAN 

266 

374 

492 

613 

867 

1140 

1440 

STD. DEV. 

43 

23 

29 

23 

29 

41 

95 

059 

220 

340 

466 

575 

825 

1070 

1350 


259 

350 

460 

575 

820 

1060 

1350 


300 

440 

570 

720 

973 

1250 

1520 

MEAN 

260 

377 

499 

623 

873 

1130 

1410 

STD. DEV. 

40 

55 

62 

84 

87 

107 

98 

060 

267 

415 

545 

675 

930 

1200 

1500 


275 

410 

550 

680 

950 

1230 

1570 


247 

350 

465 

580 

840 

1090 

1360 

MEAN 

263 

392 

520 

645 

907 

1170 

1460 

STD. DEV. 

14 

36 

48 

56 

59 

74 

84 

061 

256 

400 

540 

690 

955 

1240 

1540 


252 

325 

440 

550 

800 

1040 

1310 


300 

440 

575 

715 

1000 

1290 

1620 

MEAN 

269 

388 

518 

652 

918 

1190 

1490 

STD. DEV. 

27 

58 

70 

89 

105 

132 

161 

062 

278 

370 

482 

624 

860 

1120 

1410 


250 

343 

465 

600 

850 

1120 

1420 


239 

356 

480 

624 

890 

1180 

1480 

MEAN 

256 

356 

476 

616 

867 

1140 

1440 

STD. DEV. 

20 

14 

9 

14 

21 

35 

38 

GRAND MEAN 

280 

396 

521 

649 

904 

1170 

1460 

STD. DEV. 

48 

57 

65 

70 

73 

79 

9 
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TABLE XIV 

DYNAMIC SHEAR MODULU G', VERSUS FREQUENCY AT 40°F 


FREQUENCY (Hz) 


MIX 

(9910- ) 

5 

10 

15 

20 

30 

40 

50 

053 

1480 

1770 

2100 

2440 

3100 

3800 

4820 


1530 

1940 

2300 

2620 

3260 

3900 

4850 


1550 

1920 

2290 

2540 

3340 

4060 

4920 

MEAN 

1520 

1880 

2230 

2570 

3230 

3920 

4860 

STD. DEV. 

36 

93 

113 

110 

122 

131 

51 

057 

1540 

1880 

2220 

2560 

3260 

4000 

5040 


1480 

1790 

2130 

2450 

3180 

3880 

5000 


1560 

1930 

2270 

2590 

3280 

4060 

5020 

MEAN 

1530 

1870 

2210 

2530 

3240 

3980 

5020 

STD. DEV. 

42 

71 

71 

74 

53 

92 

20 

059 

1410 

1760 

2100 

2420 

3110 

3860 

4800 


1580 

1960 

2300 

2640 

3400 

4140 

5160 


1540 

1920 

2280 

2640 

3410 

4180 

5120 

MEAN 

1510 

1880 

2230 

2570 

3310 

4060 

5030 

STD. DEV. 

89 

106 

110 

127 

170 

174 

197 

060 

1740 

2110 

2480 

2850 

3580 

4280 

5120 


1670 

2050 

2430 

2820 

3540 

4280 

5240 


1400 

1720 

2070 

2420 

3120 

3930 

4860 

MEAN 

1600 

1960 

2330 

2700 

3410 

4160 

5070 

STD. DEV. 

180 

210 

224 

240 

255 

202 

1*4 

061 

1620 

2000 

2380 

2740 

3510 

4280 

5250 


1660 

2080 

2470 

2830 

2570 

4460 

5400 


1600 

1990 

2380 

2780 

3420 

4150 

5700 

MEAN 

1630 

2020 

2410 

2780 

3500 

4300 

5250 

STD. DEV. 

31 

49 

52 

45 

75 

156 

150 

062 

1590 

1940 

2300 

2700 

3530 

4370 

5300 


1610 

1990 

2340 

2710 

3100 

4090 

5000 


1470 

1810 

2160 

2500 

3250 

4400 

5030 

MEAN 

1560 

1910 

2270 

2640 

3390 

4190 

5110 

STD. DEV. 

76 

93 

95 

118 

140 

159 

165 

GRAND MEAN 

1560 

1920 

2280 

2630 

3350 

4100 

5060 

STD. DEV. 

89 

113 

127 

144 

161 

186 

171 
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TABLE XV 

DYNAMIC SHEAR MODULUS, G\ VERSUS FREQUENCY AT 40°F 

FREQUENCY (Hz) 


MIX 


(9910-) 

5 

10 

15 

20 

30 

40 

50 

053 

626 

832 

1030 

1240 

1640 

2040 

2440 


633 

870 

1070 

1280 

1690 

2100 

2520 


640 

890 

1100 

1310 

1710 

2100 

2490 

MEAN 

633 

864 

1070 

1280 

1680 

2080 

2480 

STD. DEV. 

7 

§ 

29 

35 

36 

36 

35 

40 

057 

606 

820 

1030 

1240 

1660 

2080 

2570 


601 

850 

1050 

1280 

1720 

2160 

2660 


730 

915 

1120 

1330 

176C 

2170 

2720 

MEAN 

646 

862 

1070 

1280 

1710 

2140 

2650 

STD. DEV. 

73 

49 

47 

45 

50 

49 

75 

059 

565 

796 

1000 

1220 

1680 

2100 

2640 


672 

910 

1120 

1340 

1780 

2200 

2680 


630 

913 

1140 

1380 

1820 

2250 

2730 

MEAN 

622 

873 

1090 

1320 

1760 

2180 

2680 

STD. DEV. 

54 

67 

76 

87 

72 

76 

45 

060 

685 

950 

1170 

1400 

1810 

2220 

2650 


660 

900 

1040 

1380 

1320 

2260 

2810 


577 

792 

1020 

1220 

1680 

2120 

2660 

MEAN 

641 

881 

1110 

1330 

1770 

2200 

2710 

STD. DEV. 

56 

81 

79 

99 

78 

72 

90 

061 

637 

849 

1070 

1290 

1730 

2160 

2670 


652 

850 

1050 

1280 

1720 

2170 

2720 


610 

800 

1000 

1210 

1610 

2020 

C"30 

MEAN 

633 

833 

1040 

1260 

1690 

2120 

2o-3 

STD. DEV. 

21 

29 

36 

44 

67 

84 

h5 

062 

625 

885 

1110 

1340 

1800 

2240 

2250 


659 

930 

1150 

1340 

1780 

2060 

2400 


609 

810 

1020 

1240 

1650 

2070 

2550 

MEAN 

631 

875 

1090 

1310 

1740 

2120 

2570 

STD. DEV. 

26 

61 

67 

58 

81 

101 

176 

GRAND MEAN 

634 

865 

1080 

1730 

1730 

2140 

2630 

STD. DEV. 

39 

50 

55 

66 

66 

74 

111 
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3.3 ROCKWELL INTERNATIONAL TESTS 

The Space Division of Rockwell International has conducted 
dynamic shear modulus tests in support of the Quarter Scale GVT, 
program [8,9]. The propellant, which is inert, is designated UTI-610A 
and was supplied by CSD. 

The specimen dimensions were 4 inch by 4 inch by 1 inch thick. The 
four inch square surfaces were bonded to aluminum plates. One of the 
plates was clamped to a rigid, fixed support and the other plate was 
subjected to -.ntrolled sinusoidal forces. Tests were conducted at 
67°F from 10 to 200 Hz with force amplitudes of +1 , + 2, +4 and +8 
pounds (i.e., average shear stresses of 0.0625, 0.125, 0.250 and 0.500 
psi). Test data was obtained from accelerometers used to measure the 
motion of the specimen along all three axes and rotation about the three 
axes. Test results are presented in Table XVI as a function of 
frequency and load and in Table XVII as u * net ion of frequency and 
dynamic shear strain level. The results show significant nonlinearities 
at the strain levels tested. 
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TABLE XVI 

UTI-610A DYNAMIC SHEAR MODULUS VERSUS FREQUENCY AT 67°F 


'Sinusoidal 

NForce 

+ 

1 lb 

+ 2 

lb 

+ T 

lb 

+ 8 

lb 

Frequency 

6* 

G" 

G‘ 

G" 

G 1 

G" 

G' 

G" 

(Hz) 










10 


2280 

805 

2070 

800 

2070 

800 

1900 

790 

20 


2980 

1100 

2730 

1000 

2600 

980 

2330 

905 

40 


3930 

1540 

3570 

1410 

3330 

1370 

2920 

1250 

70 


4800 

2000 

4310 

1830 

4010 

1750 

3470 

1590 

100 


5360 

2300 

4840 

2120 

4450 

2030 

3850 

1830 

150 


6040 

2700 

5440 

2500 

4990 

2380 

4290 

2140 

200 


6500 

3000 

5890 

2740 

5350 

2640 

4600 

2370 


TABLE XVII 

UTI-610A DYNAMIC SHEAR MODULUS, 6', VERSUS 
DYNAMIC SHEAR STRAIN AT 67°F 


frequency 







Strain'S. 

(%) 

20 

40 

70 

100 

150 

0.001 

3000 

4000 

4800 

5450 

5950 

0.002 

2950 

3850 

4600 

5200 

5550 

0.005 

2750 

3550 

4100 

4600 

4850 

0.01 

2550 

3150 

3650 

4000 

4170 

0.02 

2250 

2650 

3000 

3300 1 

1 3450 
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IV. DYNAMIC CHARACTERIZATION OF TP-H1I48 PROPELLANT 


Static and dynamic mechanical characterization tests were conducted 
by CSD under sub-contract to W. L. Hufferd & Associates. Dynamic shear, 
constant rate to failure and stress relaxation properties were measured 
on six batches of the Space Shuttle SRM solid propellant (-6, -7, -8, -9, 
-9970096, -9970115) supplied to CSD by Thiokol/Wasatch Division. The 
final report for this work [18] is included as Appendix A to this 
report. The dynamic shear modulus test results are summarized in the 
following sub-sections. 

4.1 TEST EQUIPMENT AND TEST PROCEDURE 

Small strain dynamic shear properties of solid propellants and 
other low modulus materials are routinely measured at CSD using a 
piezoelectric transducer device. Illustrations of the dynamic test 
device and the associated electrical equipment are presented in Figures 
12 through 14. The device employs a stack of piezoelectric crystals as 
tie dyiamic driver source. This dynamic strain is transmitted through 
t.-»e so’ id propellant sample as shear and the output load is measured 
jiing single piezoelectric crystal which functions as a load cell. 

The dynamic shear test device provides a direct measurement of 
dynamic properties as a function of frequency and test temperature. 
Equipment is calibrated using stainless steel rings inserted between the 
piezoelectric driver stack (Figure 13) and the readout crystal to 
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calibrate the test device ''ver a wide frequency range. During propellant 
tests the output signal and the phase angle are measured to determine 
the dynamic properties of the propellant over the desired frequency 
and test temperature ranges. 

Both the dynamic strain amplitude and the static compressive strain 
levels can be varied to determine the modulus sensitivity to various 
loading conditions. 

4.2 TP-H1148 PROPELLANT DATA 

Dynamic shear modulus behavior was measured for each of the six 
batches of Shuttle SRM propellant. Production batch No. TP-H1148-9 was 
evaluated at three static and three dynamic strain levels. Test data 
for one strain level were obtained for each batch of propellant at 
90°, 72° and 40° F. 

Dynamic modulus curves are presented in Figures 15, 16, and 17 for 
batch No. TP-H1 1 48-9 at 0.001* dynamic strain and static strain levels 
of 3%, 6% and 12%. A comparison of G' curves for each of these tests 
is presented in Figure 18. Static strain sensitivity at 100 Hz for 
batch No. TP-H1 148-9 illustrates the dynamic modulus increase with 
compressive static strain levels: 


Compressive Strain, % G* , psi 

3 2,070 

5.8 2,950 

12 6,250 


where temperature * 72° F, and dynamic strain = 0.001%. 



Dynamic shear modulus (6’ and G"). psi 
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Figure 17, Master Dynamic Shear Modulus for Batch Mo. TP-H1 148-9 at 12> 


Static and 0.00 If Dynamic Strains 
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Figure 18- Comparison of Static Strain Level* fcr Batch Mo. TP-H1 148-9 
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Test data in Figures 15, 16 and 17 shows that 6' actually drops 
lower than G" at the 90° F test temperature for the 3% and 6% static 
strains. The soft propellant at 90 C F is absorbirg a larger percentage 
of the energy input to the sample. 

Empirical shift factors for the three temperatures worked well at 
all three strain levels (Figure 19). 

Dynamic shear moduli as a function of dynamic strain level was 
measured and is presented in Figures 20 through 23. The real part of 
the dynamic modulus increases with deceasing strain level but the total 
range in modulus is less than a 50% variation. This is a much smaller 
amplitude change than experienced with static strain variations. 

Dynamic modulus test results for the other batches of TP-H1148 
propellant are presented in Figures 24 through 28. These batches were 
very similar to the major batch (-9). The real part of the dynamic 
modulus " / »r each of the six batches is presented in Figure 29. 

The real part of the dynamic modulus for each batch at 100 Hz is 
given below: 

Batch No. 6* , psi at 100 Hz 


9970115 

3,100 

-9 

2,120 

-7 

2,120 

-8 

2,100 

9970096 

1,900 

•6 

* 

CO 

o 

o 


3atch No. 9970115 has the highest dynamic modulus with the others being 
close together; batch Mo. TP-H1148-6 is the lowest. 
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Figure 20 . Master Dynamic Shear Modulus for Batch No. TP-H1143-9 at 3% 
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Figure 21. Comparison of Shear Modulus for Batch No. TP-H1148-9 at 2.9X 
Static and Different Dynamic Strain Levels 
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Figure 22. Comparison of Shear Modulus for Batch No. TP -HI 148-9 at 5.81 
Static and Different Dynamic Strain Levels 
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Figure 23. Coag>arlson of Shear Modulus for Batch Mo. TP-HL 148-9 at 121 
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Frequency (CuAy), Hi 

Figure 24. Master Dynamic Shear Modulus for Batch No. TP-H1148-6 
at 31 Static and 0.001% Dynamic Strains 



i gure 25- Master Dynamic Shear Modulus for Batch No. TP-H1148-7 
at 3% Static and 0.001% Dynamic Strains 
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Figure 26. Master Dynamic Shear Modulus for Batch No. TP-H 1148-8 
at 3X Static and 0.G01X Dynamic Strains 
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Figure 27. Master Dynamic Shear Modulus for Batch No. TP-H1148-9970096 
at 31 Static and 0.001% Dynamic Strains 
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The real and imaginary part of the dynamic shear modulus values 
at 50 Hz are given in Table XVIII for batch -9 and Table XIX for the 
other five batches. Individual values are listed along with mean 
values. 

4.3 DISCUSSION OF TEST RESULTS 

The test data obtained show very small batch-to-batch propellant 
variability and excellent sample-to-sample reproducibility. Dynamic 
modulus values are sensitive to the static compressive strain level as 
well as dynamic strain. Some fuel pockets and voids were observed in 
the small strain dynamic shear test specimens which may explain the 
larger dat< ariability than observed with the larger test specimens 
used for the other test modes. 

The test strain levels covered 0.001% to 3% and, correspondingly, 
a large modulus range. This strain range should be applicable to the 
Shuttle SRM dynamic loading conditions. Specific modulus values for 
dynamic analyses can be selected from the available data by defining 
the propellant strain range and then used with available NASTRAN 
computer analysis results. 

The compressive strain condition of the dynamic shear sample 
resembles the SRM motor loading during ignition as the propellant grain 
will be in compression from the internal pressure and a shear load will 
be superposed due to gravity, loading and acceleration. This particular 
dynamic shear test [19,20] is more appropriate for determining propel- 
lant dynamic properties than the Gottenberg disk and other dynamic 
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tests inasmuch as large strain gradients exist through the test 
sample which are not representative of the Space Shuttle SRM 
conditions. 
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V. SRM PROPELLANT DYNAMIC RESPONSE MODEL 


When considering the stress-strain relation of an elastic material, 
it is evident that for a particular value of stress there is associated 
a particular value of strain, and regardless of the length of time that 
the stress acts on the body, or what path was followed in applying it, 
the strain remains constant. In viscoelastic materials, on the other 
hand, when a stress is applied to the body, the strain state depends 
upon the manner in which the stress is applied; that is, whether the 
load is applied rapidly or slowly. Thus, the history of loading, as 
well as the magnitude of the load must be considered in describing the 
response of a viscoelastic material. In addition, a viscoelastic body 
will not maintain a constant deformation under a constant stress, 
regardless of the loading pattern; rather, it will deform or creep 
with time. Also, if such a body is constrained at constant deformation, 
the stress necessary to hold it constrained gradually diminishes or 
relaxes with time. 

The stresses and strains at a point in a viscoelastic body may 
thus be expected to vary with time, or the frequency of loading. And 
consequently, also with temperature, as will become evident. 

The following paragraphs discuss procedures for representing the 
viscoelastic response properties of solid propellants in general and 
TP-H1148 in particular. Additional details and justification of the 
applicability of the procedures to solid propellants are presented in 
References 12 to 15. The development herein is in terms of the modulus 
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of the material. Completely analogous results can be obtained employ- 
ing the creep compliance representation of the response. 

5.1 IS OTHERMAL RESPONSE 

A convenient representation for the uniaxial (or shear) stress, 
a, in a linear viscoelastic material is the so-called integral or 
relaxation formulation which, for the isothermal case, has the form 

« (t) •/ E rel (t-T) Mil dT 0) 

.00 

where E re -| ( t-x ) is the relaxation modulus in tension (or shear) and 
e(t) is the imposed strain history. 

The relaxation modulus is defined as the stress decay associated 
with a step input strain e Q . With e * c q H(t), Equation (1) becomes 

° (t) " f.Jre 1 (t ' T) e o <(t) dT (2) 

or 

2llL . E , t) ... 

e Q t rel ( 3 ) 

where H(t) is the Heaviside unit step function and 6(x) is n* Dirac 
delta function. 

For a constant strain rate test, with strain rate R (e = Rt), 
Equation ( 1 ) yields 

• ott) ■ R f 0 \el <*-*> ^ 


(4) 
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Differentiating Equation (4) with respect to time, t, gives 


RE r + 


■ R E rel (0) - R 


£ 


«rel 


■ar 


dx 


R E r,1 


or 



e ■ Rt 


(5) 


Thus, the relaxation modulus for a linear viscoelastic material can 
be deduced from differentiation of the stress-strain curve i r a constant 
rate test. 

An alternative method of measuring, and expressing, viscoelastic 
behavior is by considering steady-stc..e response to forced vibrations. 
Creep and relaxation experiments are not capable of providing complete 
information concerning the mechanical behavior of viscoelastic solids. 

In certain cases the response of a structure is sought to a loading 
for times substantially shorter than the lower limit of time of usual 
creep or relaxation experiments. 

For a linear system, botn stress and strain will vary si’-’*>oidally 
with the same frequency as the forc.ng frequency: 



a(t) = a e 


iojt 


e*(w)e 


iait 


If the input to a linearly viscoelastic material is an oscillatory 
stress. 


a<t) * ®0 e 


iwt 


( 6 ) 


ther the strain response will be an oscillation at the same frequency 
as the stress, bjt lagging behind by a phase angle 5. Thus, 

cM ■ e„ e' 1 ^- 51 (- 

where e is the strain amplitude. The phase angle, 6, is often called 
c 

the loss angle and is a function of the internel iction. It is 
.onvenient to write Equation (7) in tne form 


ei«) 


• -i<S, 

ve Q e )e 


iuit 


- t 

- ce 


(8) 


where e* is the complex strain amplitude defined by 
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If the input is an oscillatory strain. 


e(t) = e Q e 


iwt 


( 10 ) 


then the stress response will lead the strain by the phase angle 6, 

a(oj) * a e 1 ^* +<s ) * a * p iwt 

° 0 e HD 

where 

o* * a Q e 15 = a Q (cos5 + i sin6) (12) 


The complex dynamic modulus is defined as the oscillatory stress 
response to an oscillatory strain input; that is, 

E*(w) ^ = J2L !l! 

£ o 


(cos<5 + i sind) 


= E’(w) + i E M (ui) 


(13) 


where E‘ (u) is in phase with the strain and is called the storage 
modulus 

E' (w) * — C os6 (141 

e o 

For low or medium damping materials is nearly equal to the relaxation 


modulus. 
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The second term of the last line of Equation (13), E"(oj), is 
often called the loss modulus and is defined by 

E"U>) » -2- sin6 (15) 

o 

Equation v '3) may also be written in the form 

E*(u>) * |E*j a 1 5 (16) 

where 

|E*| * 
and 

6 * 

The ratio defining tan 6 is called the loss tangent or mechanical 
loss. It gives a direct measure of the energy dissipation or damping 
characteristics of the materials. For an elastic material 6=0, and 
E" = 0 giving rise to an instantaneous response. For a viscous fluid 
6 = oo and E' = 0 with the response 90 degrees out of phase. A visco- 
elastic solid has a loss tangent between these limits; i.e., 0 <6 < ®. 

A representation for the complex dynamic modulus can be obtained in 
terms of the stress relaxation modulus by substituting Equation (10) 
into Equation (2). Then, 

o(t) • s f “ f * E re1 (t-T) e iuT dT < 19 > 
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Performing a change of variables by letting u * t - x, Equation (19) 
becomes 


o(t) 



(u) e“ la)U du 


( 20 ) 


or 


E*(u>) * 


E o e 


iwt 


lb) f 

J n 


*rel 


(u) e 


-IGJU 


du 


( 21 ) 


Recalling the definition of the Laplace Transform, Equation (21) is 
seen to be 




usiw 


( 22 ) 


Ihus, if the stress relaxation modulus is known, the dynamic 
modulus can be obtained by taking the Laplace transform, setting the 
transform variable equal to (iu>) and multiplying the result by (iw). 
This observation is particularly useful if an exponential series (Prony 
Series or Dirichlet Series) is used to represent the relaxation modulus 
in this case, 

n 

E rel * E e + H E k ex P (“t/x. ) (23) 

k-1 * 


where E g is the equilibrium relaxation modulus, and E^ and are 
constants chosen to fit the experimental data. E' and E" are readily 
determined from Equation (23) to be 
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and 


E 1 (u) 


n 



E lc 

1 + («*»T k ) 2 


E" (u») » 



E k (h)T k ) 

1 ♦ («T k ) 2 


(24) 


(25) 


However, this conversion does not work well for solid propellants which 
usually exhibit nonlinear modulus variations [7,16,17]. 

Inasmuch as the conversion of static to dynamic material proper- 
ties does not work well for solid propellants and also due to the fact 
that most propellant data, including the shuttle propellant, exhibit 
log-linear behavior over most of the time or frequency range of 
interest, it is often more convenient to represent the rela>ation and 
dynamic moduli as power laws. 


E(t) = E q t" n 
E' («) * E 0 w n 

E"(u>) = EJ u> n " (26) 


where the exponent represents the slope of the modulus versus log time 
or log frequency plot. Equation (26) describes the dynamic response 
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of the space shuttle propellant particularly well since we are 
concerned only with the response over a narrow temperature and 
frequency range. 

5.2 TIME-TEMPERATURE SUPERPOSITION 

The analyses of the previous section apply to isothermal conditions, 
i.e., constant temperature conditions. In practice, it is impossible to 
obtain either relaxation data or dynamic data over the time scales and 
frequency ranges of practical interest in a solid rocket motor. In 
order to obtain data in these ranges the concept of time-lemperature 
superposition is routinely employed. (See References 12 to 15 for a 
general discussion.) 

It has been widely found that temperature has the effect of expand- 
ing or contracting the time scale of response viscoelastic materials and 
that an equivalence tetween time or frequency and temperature exists. 

Time and frequency are roughly the inverse of one another so that short 
time or high frequency response at one temperature corresponds to 
longer time or lower frequency response at a lower temperature. The 
converse holds true at higher temperatures. Thus, by obtaining the 
propellant response at several temperatures, a time-temperature shift 
function, a-p relating the equivalence of time or frequency and tempera- 
ture can be experimentally determined by horizontally "shifting" the 
test data so that it superimposes to form a single curve at some given 
reference temperature, usually 21 to 25°C. The resulting curve is known 
as the "master" response curve and is expressed in terms of temperature- 
reduced time, t/a-p in the case of the master relaxation modulus curve. 



-75- 


or in verms of temperature-reduced frequency, a>aj, in the case of the 
master dynamic moduli curves. 

Thus, the previous analyses remain valid in terms of the master 
relaxation and dynamic moduli if time, t, and frequency, u>, in the 
previous expressions are replaced by temperature-reduced time, t/a T , and 
temperature-reduced frequency, wa-p respectively. 

A common analytical representation of the shift factor for amorphous 
polymers is known as the WLF Equation, which is based on free volume 
considerations in the neighborhood of the glass transition temperature, 

Tg [2U : 

-8.86(T-T ) 

log a T = 5 (26) 

1 101.6 +T-T $ 

where T is measured in degrees Kelvin. 

A somewhat better representation of the temperature shift factor 
for solid propellant is obtained using the modified power law represen- 
tation [14] 



where is the reference temperature, and is a characteristic 
temperature. The representation given by Equation (27) is more conven- 
ient for analyzing propellant test data than the WLF equation. 
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VI. ANALYSIS OF VARIANCE MODEL 

The problem of variability in solid propellant mechanical properties 
has been known for a long time. Bills and Svob [22] discussed the 
"large statistical variation" in a composite propellant mechanical 
properties as i relates to establishment of mechanical properties 
criteria. Marti, Morill and Bersche [23] reported typical values of 
the standard deviation for various material properties studied during 
an aging program. Majerus [24] presented histograms of strain data 
from uniaxial and biaxial tests on propellant, and Wiegand [25] corre- 
lated maximum stress and modulus for uniaxial tensile samples tested 
under various conditions of humidity and temperature and found a strong 
positive linear correlation between stress and modulus when plotted on 
log-log paper. An application of statistics to failure criteria was 
given by Briar and Wiegand [261, and in another publication, Briar [27] 
presented application of Monte Carlo simulation in modeling constant 
strain failures in solid Dropellants. More recently statistical tech- 
niques have beer, applied to the structural analysis of a cast double-base 
solid propellant grain [23] and the interpretation of elevated temperature 
data in connection with propellant service life predictions [29]. 

Many source* . f variation operate together to produce the final 
variability in solid propellant mechanical properties. Variability in 
propellant ingredients, curative ratio, temperature and duration of cure, 
test sample machining, storage, handling, testing and data reduction are 
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all responsible for the variability in the final propellant test data. 

This variability may be classified in two major components: 

(1) variability due to inhomogeneities in the 
propellant (propellant variance). 

(2) variability from all other sources including 
testing, machining, storage, handling, data 
reduction, etc. (sampling variance). 

The first component, propellant variance, should be used in a statistical 

stress analysis since it represents Ine actual variability of the 

physical property of interest. The second component does not characterize 

the propellant, but occurs because the sampling procedure is imperfect. 

The objective of the statistical analyses of the SRM propellant test 

data is to determine the form of the probability distributions, the 

applicable parameters of the probability distributions and the extent, 

if any, of correlations between parameters. 

Previous experience with solid propellants indicates that most data 

is normally distributed and thus the distribution may be characterized 

by mean and variance. The applicability of a normal distribution may 

2 

be investigated from cumulative frequency diagrams using the x -goodness- 
of-fit test, and by investigating the asymmetry or skewness and the 
kurtosis of the assumed normal distribution [30]. If the form of the 
distribution is not normal, then other closed form probability distribu- 
tions may be considered; e.g., log normal, exponential, gamma, etc. If 
an exact closed form probability distribution cannot be obtained, then 
an approximation to the distribution may be generated. 
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The following subsections present the analysis of variance for the 
two ma.ior components of propellant variability and a probabilistic 
analysis of the SRM propellant dynamic response model. These analyses 
are applied to the SRM shuttle propellant in the next section of this 
report. 

6.1 ANALYSIS OF VARIANCE 

An analysis of variance is a technique by which the variations 
associated with certain components of the variance may be isolated and 
estimated [30,31 ]. 

When two or more independent sources of variation operate, the 
resulting variance is the sum of the separate variances. The procedure 
followed here is for two independent variances: 

2 

(1) variance between castings denoted by 

(2) variance within a given batch denoted by a Q 

The total variation is obtained by simple addition of the two independent 
components: 


a 2 - a 2 + a 2 (28) 

To aetermine the variance exactly requires an infinite number of 

observations; thus, in practice the variance may only be estimated. 

If, in the statistical model, it is assumed that test samples from 

each batch can be regarded as random elements from a normally distributed 

population of test samples whose variation is completely specified by 
2 

the variance o Q , then each casting may be regarded as a random element 
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from a normally distributed population of batches whose variation is 

2 

completely specified by the variance o-j. The purpose in analyzing the 

2 2 

data is to estimate a Q and a-j. 

Suppose that there are k batches and n propellant samples tested 
from each batch; then the total number of samples is N = kn. Denote the 

observations on the i th batch by (X^, X^, X^ » X^) with mean X.. 

An estimate of the within-batch variance is then given by: 

<29) 

where the arrow means "is estimated by". A similar expression can be 

obtained for each sample, and on the reasonable assumption that the 

2 

variance does not vary from batch-to-batch, a better estimate of a Q may 
be obtained by averaging the variance e c ti...«ted from each batch: 



If the means of batches a> i denoted by (X-j, X^, - - -, X^) and 
the grand mean by then an estimate of the variance of the mean is 




( 31 ) 


The variance of the mean of n tests on a given batch due to within-batch 
2 

variation is a . The k batch means will differ also because of the 

2 

batch-to-batch component of variations and since the within and 
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between sources of error are independent, tneir variances are additive 

2 2 

and the variance of each sample mean is + cr/n, or 

a l + a o /n " FT .£ (X i " X)2 (32) 

To facilitate the Analysis of Variance and to help organize the 

calculations, the format in Table XX is followed. Table XX indicates 

that the Sum of Squares and Degress of Freedom are additive so that 

S-j + S Q = S and (k-1) + k(n-l) = n(k-l). Solving the expressions given 

2 2 

in Table XX results in estimates of n and a,. 

o 1 

a 2 + M (33) 

0 0 

A "■ ( M l " M 0 )/n ( 34 > 

TABLE XX 

ANALYSIS OF VARIANCE TABLE FOR TWO SOURCES OF VARIATION 


Source of 
Variation 

Sum of Squares 



Degrees 

of 

Freedom 

Mean 

Square 

Estimated by 
Mean Square 

Between 

lots 

LO 

II 

04 

• X 
• 

1 X*" 

c 

k-1 

St 

m i " FT 

2 . 2 
°0 * n °1 

Within 

lots 


k(n-l) 

M - S ° 

M o k(n-l) 

o 2 

0 

Total 

t £(x 13 -x) 2 - s 

i=l j*l 1J 

nk-1 
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The previous procedure may be applied when there are more than two 
sources of variation. For example, it is possible to introduce variations 
in propellant ingredients such as oxidizer lot-to-lot variations, polymer 
lot-to-lot variability, curative lot-to-lot variability, etc. The 
analysis of variance is not difficult; however, a very substantial test 
program is required for such an analysis. Such data does not exist; 
therefore, we have chosen to investigate only the influence of batch-to- 
batch variability and within-batch variability on the SRM propellant 
dynamic response. 

6.2 CORRELATION ANALYSIS 

A correlation analysis of the relationship between two variables 
begins with an attempt to discover the approximate form of the relation- 
ship by graphing the data in a so-called "scatter diagram"; a plot of the 
points representing corresponding pairs of measurements of the two 
variables being studied. The most common relationship is a linear 
relationship where the data tend to form a straight line. 

The usual measure of the degree to which the variables are linearly 
related is the correlation coefficient r. For two variables, X and Y, 
the correlation coefficient is defined by; 

r " no^Oy (X i' X)( Y^ (35) 

The absolute magnitude of r indicates the degree of correlation between 
the two variables. When all of the points lie on a stio^ght line 
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(perfect linear correlation) the correlation coefficient is +1. T^e 
correlation coefficient is zero when there is no linear correlation 
between the two variables; however, there may be a nonlinear correlation. 
A positive correlation coefficient inoicates that an increase in one 
variable tends to be accompanied by an increase in the second variable, 
and a r.wgative correlation indicates an inverse relationship. 

It is also useful to consider the correiation coefficient with 
record to a linear regression model. In using a lir ■* regression 
model with one dependent variable (X). it is assumed that the inde- 
pendent variable (Y) is normally distributed with constant variance 

2 2 
o X y at each value of the dependent variable. Thus, Oy X is an 

estimate of the variation of Y about the regression line Y on X. The 

square of the correlation coefficient indicates the relative portion 

of the variation in Y which can be explained by the dependence on X, 


2 

a YX 


( 24 ) ( 1-44 


(36) 


For a multivariate correlation analysis of several variables, 
such as the situation with the SRM propellant dynamic response model, 
the covariance, Eouation (36) becomes a matrix of the form 
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Confidence limits for the estimates for the means, standard 

deviations and correlation coefficients can be estimated using the 
2 

X goodness-of-fit test or the students t-distribution. 

6.3 PROBABILITY ANALYSIS OF SRM PROPELLANT DYNAMIC RESPONSE MODEL 

The SRM propellant dynamic response model is a power-law relation- 
ship of the form 


E = E 0 (wa T ) n 


(38) 


for either the real (6' or E*) or imaginary (G H or E") parts of the 
complex tensile (E*) or shear (G*) moduli, where u is the circular 
frequency and a T is the time- temperature shift factor defined by 


a 


T 



(39) 


and T r is the reference temperature to which data at other temperatures 
are shifted and T and m are experimentally determined parameters. 

a 

The constant E Q in Equation (38) represents the value of the dynamic 
modulus at a temperature-reduced frequency of ua T = 1, and n is the slope 
of the curve on a log E versus log wa T plot. 

A computer code, developed in [4] performs a least-squares 
curve-fit of the master dynamic modulus versus temperature-reduced 
frequency test data to determine the parameters E q and n, or 
predicts dynamic modulus as a function of temperature and frequency, 
given E Q and n. The computer code is applicable to the real and 
imaginary parts of the t.nsile or shear modulus. 
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In order to account for propellant variability it is desirable 
to incorporate probability concepts into the SRM propellant dynamic 
response model. There are basically four approaches for accomplishing 
this objective [32]: 

(1) Monte Carlo Simulation 

(2) Maximum Likelihood Methods 

(3) Moment Generating Functions 

(4) First-Order, Second -Moment Theory 

Monte Carlo simulation techniques consist of evaluating Eq. (38) 
using randomly selected values of the input parameters, E Q and n. In 
order to use the method and select representative values of E Q and n 
at random it is necessary to first have a statistical description of 
each parameter defining the manner in which the parameter can be expected 
to vary. 

Maximum likelihood methods deal primarily with finding a point 
estimate of the variance nased on multiple sets of values for a function 
of more than one variable. For functions of one variable, this method 
of estimation is efficient and with large sample sizes is very consis- 
tent. If the functional relationships are known, the likelihood methods 
may be extended to functional combinations. 

The variance of a function h(x^,x 2 ... xj of n independent 
random variables may be computed by means of moment generating functions. 
Let 


z * hU^x^ ... , x n ) 


(40) 
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_ 2 

then the variance of z is given by the expected value of (z-z) , 
E[(z-z) 2 ], where 

E[(z-i) 2 ] ■ f “<z-z) 2 f(z)dz (41) 

^-00 


and f(z) is the probability density function of z. 
Therefore, 


E[(z-z) 2 ] * 


r f 


- 2 

[h(x-| ,x^ , • • • ,x )-z] 


* f(x^ .... ,x n )dx-j . . .dx r 


(42) 


Now a moment generating function for z is given by 


MGF(z-z) 


• f° e t(2 ‘ J) f ( z ) dz 
J *00 


(43) 


Differentiating twice with respect, to t and setting t = 0, the variance 
of z is obtained: 


var[z] 


d MFG(z-z) 

— — L 

dr 



(z-z) 2 f (z)dz 


(44) 


- * , t*0 

Using suitable transformations one can then transform to integrals 
involving x^ .Xg,. . . ,x n - 

First-order, second-moment theory probably offers the most 
advantage because of its relative ease of application and accuracy. 
Cornell [33] in attempting to encourage wider usage of probability 
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concepts by non-specialists, proposed a so-called "first-order second- 
moment" reliability theory. 

The term "first-order" indicates the linearization of 
nonlinear equations and the term "second-moment" refers to 
the fact that only the first two moments of the involved random 
variables are needed. For example, some random variables can be 
characterized by their mean X and variance a x (or equivalently their 
standard deviation, o , or coefficient of variation V * a /X). The 

A A A 

stochastic dependence between a random variable X and a random variable 
Y is measured by their covariance Cov [X,Y] (or equivalently by their 
correlation coefficient, r X y = Cov[X,Y]/a x o y ). 

The application of first-order, second moment theory is illus- 
trated in the following example. Let Y be a continuous function of X, 
i.e., Y = g(X), where X and Y are random variables. Denote the 
mean of X by X or E[X]. Taking a Taylor series expansion of g(X) 
about X,g(X) may then be expressed as: 


g(x) • g (x ) * (x-x) 4am 




(45) 


Keeping only the first two terms in the expansion and taking the 
expectation of both sides. 


E[Y] = E[g(X)] = g ( E[X] ) 


(46) 


since E(X-X) * 0. Similarly, keeping the same terms and finding the 
variance of both sides yields the approximation 
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I *\ 2 

X xl 


Since g(X) is a constant. 


(47) 


Var[g(X)] * 0 (48) 

and 



» 

T An { 

* 

y \ I T - TAni 

r y \| 


Var 

ra 

^ l (x ‘ x) * 

] Var[X] 

(49) 


L 

m 

1 X-X J J ^ 

“ 1 

1 X 
N 

X 



Using this same principle for a variable, say Z, which is a 
continuous function of several variables, X^, X 2 , . .., X n (whose means, 
variances, and covariances are known), the first order approximation 
of the mean and variance of Z can be obtained by expanding the function 
g(X 1 ,X 2 ,... ,X n ) in a Taylor series about the mean of the X's. Keeping 
only the linear terms in the deviation from the means and applying 
the same linear transformation rules for the mean and variance, yields 


z * g(x r x 2 x n ) 


and 


Var[Z] 



+ 




(50) 


(51) 
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in which the derivatives are evaluated at the mean values of X.. and 

V 

If the random variables X^ and Xj are independent, then Eq. (51) 
simplifies to 


Var[Z] 



(52) 


Treating E„, n, T and m as random variables in Eqs. (38) and 

0 a 

(39) and ui, Tp and T as deterministic, then 


E = E(E 0 ,n,T a ,n) 


and 


1L 

3E, 


E=E 


/ “ \ n E 

= {vAj) 3 — 

E o 


3E 

3n 


_ - M 

E*E W *T 


IL 

3T, 


Ie*e 


3E aa T 
3a T 3T a 


E*E 


- (T-Tp) 
n m E — x- 


dE 

3m 


E*E 


aE_!!i 

aaj am 


. 3 n m E 


E*E 


li 

V^a 


(53) 

(54) 

(55) 

(56) 

(57) 


Applying the previous developments, i.e., 
follows that 


Eqs. (50) and (51 ) , it 
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- a nm n m 

u>aj 
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r /t t \2 fnT 
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If the random variables E_, n, T and m are not 

0 d 


(58) 

(59) 



correlated 


then Eq. (60) reduces to 
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Var[E] = o 2 = (nE) 2 


.0 (T-Tr) 

+ m“ - 


) o 

w 


n 

77? 


- — r a T + in 

-r \ ^ 


( 3 t ) 


2 l J -\ 


<VV * 


<VV 


,1 

m « 


( 61 ) 


A distinct advantage of this approximate closed form proba- 
bilistic approach is that the results are independent of any assumptions 
regarding the form of the underlying probability distributions. 
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VII. ANALYSIS OF SRM AND RELATED PROPELLANT DATA 


Statistical analyses of UTI-610 inert propellant test data and 
SRM propellant test data generated by Thiokol and CSD have been 
analyzed according to the statistical taodels of Section VI in order to 
estimate the expected variance of propellant dynamic response proper- 
ties in the Shuttle SRM. In addition, the influences of dynamic and 
static strain levels are considered and incorporated as necessary into 
the SRM propellant dynamic response model. Based on accelerated aging 
test data of Thiokol [10] an estimate of aging behavior is also made. 

7.1 STATISTICAL ANALYSIS OF PROPELLANT DATA 

The probabilistic analysis developed in Section 6.3 has been 
applied to the University of Utah test data on UTI-610 inert propellant 
in an attempt to estimate the influence of model parameters on dynamic 
response. An analysis of variance was conducted on the Thiokol and 
CSD test data to determine expected within-batch and batch-to-batch 
visibility. 

The results of statistical analysis of UTI-610 propellant are 
presented in lables XXI through XXVI. Tables XXI and XXII present the 
model parameters at -18°C, 4°C, 24°C and 49°C and static strain levels 
of 0.1, 0.5, 1.0, 2.5, and 5.0 percent. The temperature shifted 
model parameters are presented in Tables XXIII and XXIV. Substantial 
dependence on static strain level is noted; however, the coefficients 



TABLE XXI - MODEL PARAMETERS, REAL PART OF DYNAMIC TENSILE MODULUS AS A FUNCTION OF TEMPERATURE 
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TABLE XXII - MODEL PARAMETERS, IMAGINARY PART OF DYNAMIC TENSILE MODULUS AS A FUNCTION OF TEMPERATURE 

AND STRAIN 
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TABLE XXIir - TEMPERATURE-SHIFTED MODEL PARAMETERS, 
REAL PART OF DYNAMIC TENSILE MODULUS 




AS A FUNCTION OF STRAIN 


Strain 

e ; 

\ 

n 

a n 

r 

(%) 

(MN/m 2 ) 

(MN/m 2 ) 




0.1 

16.8 

1.41 

0.056 

0.0025 

0.88 

0.5 

8.76 

0.24 

0.22 

0.0032 

0.99 

1.0 

8.20 

0.26 

0.25 

0.0040 

0. .4 

2.5 

6.83 

0.31 

0.28 

0.0058 

0.97 

5.0 

5.56 

0.58 

0.28 

0.014 

0.88 


TABLE XXIV - TEMPERATURE-SHIFTED MODEL PARAMETERS, 

IMAGINARY PART OF DYNAMIC TENSILE MODULUS 
AS A FUNCTION OF STRAIN 


Strai n 

E ; 

Q 

m 

o 

n 

a n 

r 

(X) 

(MN/m 2 ) 

(MN/m 2 ) 




0.1 

3.88 

0.26 

0.36 

0.01 

0.96 

0.5 

3.43 

0.14 

0.33 

0.0072 

0.97 

1.0 

2.89 

0.13 

0.37 

C.0075 

0.97 

2.5 

2.67 

0.15 

0.38 

0.010 

0.96 

5.0 

1.72 

0.16 

0.43 

0.017 

0.91 
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TABLE XXV - MEAN AND STANDARD DEVIATION, 

REAL PART OF DYNAMIC 
TENSILE MODULUS 

STRAIN 

0 . 5 % 1 % 2 . 5 % 5 % 

E a E Q E a E 0 E a E Q E a E Q ^ 

(MN/m 2 ) (MN/m 2 ) (MN/m 2 ) (MN/m 2 ) (MN/m 2 ) (NM/m 2 ) ( NM/m 2 ) (NM/m 2 ) 


0.01 

3.13 

1.34 

2.59 

1.21 

1.91 

1.19 

1.35 

1.10 

0.1 

5.20 

2.18 

4.61 

2.10 

3.63 

2.18 

2.96 

1.80 

1.0 

8.63 

3.62 

8.20 

3.74 

6.92 

4.15 

5.63 

3.42 

3.5 

11.4 

4.77 

11.2 

5.11 

9.83 

5.90 

8.0 

4.86 

10.0 

14.3 

6.02 

14 .u 

6.65 

13.2 

7.91 

10.7 

6.52 

11.0 

14.6 

6.14 

14.9 

6.81 

13.5 

8.13 

11.0 

6.70 

35.0 

18.9 

7.92 

19.9 

9.1 

18 7 

11.2 

15.2 

9.26 

100.0 

23.8 

9.98 

25 . r 

11.8 

25.1 

15.1 

20.5 

12.4 

110.0 

24.3 

10.2 

26.6 

12.1 

25.8 

15.5 

21.0 

12.8 

350.0 

31.3 

13.2 

35.5 

16.2 

35.7 

21.4 

29.0 

17.6 

700.0 

36.5 

15.3 

42.2 

19.2 

43.3 

26.0 

35.3 

CM 

1000.0 

39.4 

! 6.5 

46.1 

19.2 

47.9 

28.7 

39.0 

23.7 
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TABLE XXVI - MEAN AND STANDARD DEVIATION, 
IMAGINARY PART OF DYNAMIC 
TENSILE MODULUS 

STRAIN 


Temperature- 

Reduced 

0.5* 

E 0- 

F 

1* 

o 

2, 

r 

.5* 

o 

F 

5* 

a 

Frequency 

(wa T ) 

(MN/m 2 )(MN/m 2 ) 

2 E 2 
(MN/nf)(MN/nf) 

E 

(MN/m 2 )'MN/m 2 ) 

2 E ‘ 
(MN/nr ) (MN/m‘ 

O.Ol 

0.74 

0.41 

O.So 

0.34 

0.47 

0.33 

0.24 

0.34 

0.1 

1.59 

0.88 

1.24 

0.73 

1.13 

0.67 

0.64 

0.77 

1.0 

3.40 

1.88 

3.91 

1.70 

2.72 

1.61 

1.71 

2.10 

3.5 

5.15 

2.84 

4.63 

2.71 

4.37 

2.59 

2.94 

3.56 

10.0 

7.28 

4.02 

6.83 

3.99 

6.52 

3.86 

4.62 

5.60 

11.0 

7.51 

4.14 

7.08 

4.13 

6.76 

4.00 

4.81 

5.83 

35.0 

11.0 

6.07 

10.9 

6.34 

10.5 

6.20 

7.91 

9.59 

100.0 

15.6 

8.58 

16.0 

9.36 

15. b 

9.30 

12.4 

25.1 

110.0 

16.1 

8.86 

16.6 

3.69 

16.2 

9,61 

12.9 

15.7 

350.0 

23.5 

13.0 

25.5 

14.9 

25.2 

14.9 

21.3 

25.8 

700.0 

29.6 

16.3 

32.9 

19.2 

32.8 

19.4 

28.1 

34.8 

1000.0 

33.3 

18.4 

37.6 

21.9 

37.5 

22.2 

33.4 

40.6 
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of variation associated with the coefficient of the power law, E q , and 
the slope, n, of the dynamic moduli curves are generally considerably 
less than 10 percent, and the coefficients of correlation are 
generally good. Even though the individual variances in the model 
parameters is small, the overall effect on the propellant dynamic 
response, as shown in Tables XXV and XXVI, is to produce expected 
variances of + 50 percent. 

An analysis of variance of variance was conducted on the Thiokol 
test data of TP-H1148 presented in Tables X through XV in Section 3.2. 
The analyses were carried out at frequencies of 10 and 50 Hz, and the 
results are summarized in Tables XXVII through XXXVIII. The value F 
given in the tables is determined by the ratio of the mean squares of 
the batch-to-batch variations and represents the F-distribution. At 
a significance level of 1 percent for the specified degrees of freedom, 
F^*^ (5,12) = 5.06, indicating that there is no significant difference 
among the means of the batches tested. 

The total standard .aviations are summarized in Table XXXIX. In 
the case of the loss modulus, G", it can only be concluded that the 
total variation due to batch-to-batch variation and within batch 
variation increase with testing frequency. In the case of G' , the 
deviation is approximately the same at 9 0 and 70°F with an increase 
at 40°F. The deviation also increases with increasing test frequency. 

The test data obtained by CSD on batches TP-H1 1 48-6, -7, -8, 

-9, 9997-0096 and 9991-0115 and summarized in Tables XVIII and XIX were 
also analyzed. The AN0VA tables are presented in Tables XL through XLV. 
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TABLE XXVII. ANOVA Table G 1 at 90°F and 10 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean Standard 
Square Deviation 


Batch-to 

Batch 

19111 

5 




Within-Batch 

40851 

! 12 

! 

3404 

58.3 


Total 

, 

59962 

17 

t 


59.5 



TABLE XXVIII. ANOVA Table G' at 90° F and 50 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean Standard 
Square Deviatior 

F 

Batch-to- 

Batch 



115178 

5 

23036 56.0 

t 

\ 

1.69 

Within-Batch 

163600 
1 

12 

! 

i 13633 , 117 

i 

! ! 


Total 

j 

278778 

l 

17 

130 

1 
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TABLE XXIX. ANOVA Table G' at 70°F and 10 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

Standard ! F 

Deviation 

j 

Batch-to- 

Batch 

11859 

5 

2372 

48.7 0.47 

\ 

1 

I 

Within-Batch 

60887 

12 

5074 

— — 1 

30.0 

Total 

72746 

17 


57.2 

TABLE XXX. 

ANOVA Table G’ at 

70°F and 50 Hz 

Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean Standard F 

Square Deviation 

Batch-to 

Batch 

75730 

5 

15150 

123 0.78 

Within-Batch 

231867 

12 

19320 

373 

Total 

307600 



129 


I 
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TABLE XXXI. ANOVA Table G' at 40°F and 10 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

Standard 

Deviation 

F 

Batch-to- 
Batch | 

55933 

5 

I 

! 

11190 

1 106 

0.84 

Within-Batch ! 

i 

i 

160067 

1 

12 

: j 

i 

13339 

\ 

26.8 

I 


i 

i 

Total 

21600 

i 

! 

17 


109 


TABLE XXXII, 

ANOVA Table G' 

at 40°F 

and 50 Hz 


Source 

of 

Variation ■ 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

J 

Mean Standard 
Squares Deviation 

. 

F 

, 

Batch-to- 1 

Batch 

i 

; 

240361 



' 

, 

1 5 

I 

. . _ 

48072 



94.0 

2.23 

! 

Within-Batch 


i 

1 12 

i 

i 

21583 

j 

i 147 


Total 

499361 

17 



174 
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TABLE XXXIII 

. ANOVA Table 6" 

at 90° F 

and 10 Hz 


Source 

of 

Variation 

Sum i 
of j 
Squares ; 

Degrees 

of 

Freedom 

Mean 

Square 

Standard 

Deviation 

F 

Batch -to 
Batch 

I 

> 

7764 j 

i 

l .. | 

5 

1553 

6.6 

1.09 

Within-Batch 

17055 

i 

12 

1422 

37.7 


Total 

24819 

17 


38.3 


TABLE XXXIV 

. ANOVA Table G" 

at 90°F 

and 50 Hz 


Source 

of 

Variation 

Sum Degrees 
of of 

Squares Freedom 

Mean 

Square 

Standard 

Deviation 

I 

F 

Batch-to- 

Batch 

59395 

5 

11879 

43.8 

1.84 

Within-Batch 

73567 

12 

6131 

78.3 



Total 


132962 17 


89.7 
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TABLE XXXV. ANOVA Table G" at 70°F and 10 Hz 


Source 

or 

Vari ation 

Sum 

of 

Squares 

! 

Degrees 
I of 
freedom 

.L- . .... _ 

Mean 

Square 

Standard 

Oeviation 

1 1 

F 

i 

i 

Batch-to- 

Batch 

j 

33802 

I 

1 5 

] 

6760 

i 

j 

. 40.4 

i 

1 

3.64 

Within-Batch 

j 

22314 

I 

12 

1860 

43.1 

| 

i 

j 

Ttoal 

1 

56116 

! 

17 


59.1 

j 

\ 

\ 






TABLE XXXVI. 

ANOVA Table G" 

at 70°F and 50 Hz 


Source 

of 

Variation 

ir 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

Standard ? 

Deviation 

( 

l 

1 

Batch-to- 

Batch 

17317 

5 

3463 

58.9 

0.36 

! 

Wi thin-Batch 

115/33 

12* 

9644 

45.4 

i 


Total 


133050 17 


74.3 



TABLE XXXVII. ANOVA Table G" at 40°F and 10 Hz 


Source 

of 

Variation 

; Sum Degrees 

S of of 

; Squares Freedom 

1 

Mean 

Square 

Standard 

Deviation 

F 1 

Batch-to- 

Batch 

4333 

5 

867 

i 29.4 

I 

0.28 


1 

1 



1 

i 


Within-Batch 

i 

' 37377 

12 

3115 

27.4 

| 

Total 

41710 

17 


40.2 

- - - - 

TABLE XXXVIII. ANOVA Table G 

" at 40°F and 50 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean Standard 
Square Deviation 

F 

Batch-to- 

Batch 

109427 

1 5 

1 

21886 

67.1 

2.61 

Within-Batch 

100533 

12 

8378 

91.5 



Total 


209961 17 


113 
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of Total Standard Deviations 


G" 
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TABLE XL. 

ANOVA Table G'.at 
0.001% to Dynamic 

90°F, 2.9% 
Strain and 

Stati c 
50 Hz 

Strain, 

Source 

of 

Variation 

Sum Degrees 

! of of 

: Squares Freedom 

Mean Standard 
Square Deviation 

■ i 

F 

Batch-to- 

Batch 

j 

; 2360544' 5 

i i 

- | . r 

> 1 

{ ! 

i 472109 ! 

1 

■ i 

93.1 

112 

Within-Batch 

t f 

i ; 

202390| 48 

J i 

j 

i 

4216 i 

! i 

i 

64.9 

i 

i 

' .... 1 


i 

Total '2562934 i 53 ; '113 


TABLE XLI . ANOVA Table G' at 70°F, 2.9% Static Strain, 
0.001% Dynamic Strain and 50 Hz 


Source 

of 

Variation 

Sum Degrees i 

of of 

1 Squares Freedom 

Mean ! 
Square I 

standard 

Jeviation 

F 

Batch-to- 

Batch 

, 7803855! 5 

« 1 

1560771 

161 

178.5 

Within-Batch 

..j | 

! t 

4721521 54 

8744 

93.5 


Total 

1 

82760061 59 

: i 

; 1 


186 
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TABLE XLII. ANOVA Table S' at 40°F, 2.9% Static Strain, 
0.001% Dynamic Strain and 50 Hz 


Source 

of 

Variation 

Sum Degrees 

of i of 

Squares Freedom 

! 

Mean S 
Square 

tandard F 

Deviation 

I 

Batch-to- 

Batch 

h ! 

* 

233157028 5 

1 

46631406 

929 840 

Within-Batch 

! | ! 

2664118 48 

I 

l 

. - l 

i 

55502 

i 

i 

236 

Total 

r j l 

235821146 53 

[ , i 1 ; 


958 


TABLE XLIII. ANOVA Table G" at 90°F, 2.9% Static Strain, 
0.001% Dynamic Strain and 50 Hz 


Source 

of 

Variation 

Sum 
i of 
Squares 

1 1 

Degrees 
of I 

Freedom j 

Mean Standard F 

Square Deviation 

1 

l. . ... .. ! 

Batch-to- 

Batch 

! 1 

' 

1591385, 

. .] 

5 ! 

318277 j 76.5 119 

! 

r ■ 1 

Within-Batch 

! 128037 

48 

! 

; i 

2667 j 51.6 

i 

i 

1 

Total 

1719423 

53 

! i 

92.3 

i 

i 

1 
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TABLE XLIV. ANOVA Table 6" at 70°F, 2.9% Static Strain, 
0.001% Dynamic Strain and 50 Hz 


Source 

of 

Variation 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

Standard 

Deviation 

I 

F 

Batch-to- 

Batch 

j 2387931 

i 

1 

1 

5 

443272 

i 

85.6 

140 

Within-Batch 

: 171510 

i 

1 

i 

( 

1 

54 

i 

[ 

i i 

i 

1 

3177 

56.4 


Total 

i 

! 2387931 

1 

! 59 

f 

■ 

j 

| 102 

1 

1 


TABLE XLV. 

ANOVA Table 6" at 40°F, 2.9% Static 

Strain, 


0.001% Dynamic Strain and 50 Hz 


i 

Source Sum , Degrees i 

of of | of 

Variation Squares | Freedom 

Mean Standard 
Square Deviation 

F 

Batch-to- 

Batch 

— - , . 

22672440 

5 

4534488 , 289 

282 

Within-Batch 

771010 

48 

j 

16063 » 127 

1 

t 

1 

Total 23443449 

53 

t 

316 j 
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Additional test data at 70°F (Tables XLVI and XLVII) were also 
analyzed in order to investigate variance as a function of static and dynamic 
strain level and frequency. The ANOVA tables for the test data in Tables 
XLVI and XLVII are presented in Tables XLVI I I through L. For the CSD test 
data we cannot conclude that the batch-to-batch means are the same. 

A summary of the total standard deviations is presented in Table LI. 
These data indicate a slight increase in the total variation between 90 and 
70°F and a significant increase between 70°F and 40°F. Tf ta presented 
in Tables XXVIII and XLVI indicate no significant depende. n dynamic 
strain level between 0.001% and 0.01% dynamic strain. A decrease in 
dynamic modulus is noted at a dynamic strain of 0.12%. 

7.2 INFLUENCE OF STATIC STRAIN LEVEL ON PROPELLANT RESPONSE 

The test data previously presented indicate considerable stiffening 
with increasing compressive strain between 3% and 12% static strains. At 
70°F the real part of the modulus approximately doubles between 3% and 6% 
and doubles again between 6% and 12%. Thus it is possible to construct a 
static strain shift factor analogous to the time-temperature shift 
factor by the relation 

0.1 5c c 

a £ * 0.65 e 5 (62) 

where a £ * 1 at 3% strain. Equation (62) applies to the real part of the 

dynamic modulus, G'. A similar analysis for the loss modulus, G", yields 

0.1 Oe 

a £ * 0.74 e s (63) 

where a £ * 1 at 3% strain and e s is again expressed in percent 


strain. 
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TABLE XIV1II. ANOVA Table G' at 70°F, 2.9% Static Strain, 
0.01% Dynamic Strain and 10 Hz 


Source 

of 

Variation 

. 

Sum 

of 

Squares 

Degrees 

of 

Freedom 

Mean 

Square 

' ' 1 

Standard i P 
Deviation 

i 

Batch- to- 
Batch 

330675 

i 3 

! 

110225 

50.3 

- 

1?.2 

Within-Batch 

326070 

; 

36 

■ ! 

9058 

95.2 


Total 

• 656745 

i 

i 

39 


108 

J 

r— — 


TABLE XLIX. ANOVA Table G 1 at 70°F, 2.9% Static Strain, 
0.01% Dynamic Strain and 100 Hz 


Source 

of 

Variation 

Sum 

° f 

Squares 

Degrees 

of 

Freedom 

Mean Standard 
Square Deviation 

F 

Batch- to - 
Batch 

1191539 

3 

! 

1 

i . 

399180 97.9 

i 

I , 

25.0 

Within-Batch 

573798 

i 

36 1 

i 

15939 

1 

126 

| 

j 

Total 

j 

1771338 

i 

39 


160 
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TABLE L. ANOVA Table 6" at 70°F, 2.9% Static Strain, 
0.01% Dynamic Strain and 10 Hz 


Source 

of 

Variation 

mm 

mm 

Degrees 

of 

Freedom 

Hear 

Square 

Standard 

Deviation 

F 

Batch-to- 

Batch 

90491 

3 

30164 

24.6 

5.97 

Within-Batch 

217599 

| 36 

1 

| 

; 6044 

! 

i 

t 

| 77.7 


Total l 

308090 

f 

1 

l ” 

i 

! 

! 

i 

\ 

| 81.5 

t 
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TABLE LI. SUMMARY OF TOTAL STANDARD DEVIATIONS AT 2.9% 

STATIC STRAIN 


Temperature 

m 

Dynamic 

Strain 

(%) 

G' 

G" 

10 Hz 

50 Hz 

100 Hz 

10 Hz 

50 Hz 

90° F 

0.001 


113 



92.3 

70°F 

0.001 


186 



102 


0.001 

108 


160 

81.5 



40°F 


O.OOl 


958 


316 
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7.3 INFLUENCE OF AGING CN PROPELLANT RESPONSE 

Solid propellants typically experience changes in chemical and 
physical properties due to normal aging during long term storage anj 
thus the SRM propellant may be expected to undergo aging. Thiokol [10] 
has conducted preliminary accelerated aging tests on TP-H1148 
propellant at 75, 110, 135 and 150°F, and interpreted the data accordinn 
to a previously developed agit.c, model [34,35]. Thiokol has shown that 
the reaction responsible for the propellant cure is essentially complete 
in one week and contributes little to the mechanical property changes 
with age. The aging rate is determined by following tht change in gel 
formation in the binder, which follows an Arrehenius rate equation. 

The initial aging data obtained on TP-H1148 obtained by Thiokol, 
although difficult to interpret quantitatively due to initial softening 
of the propellant rather than the more typical continued po c t curing, 
qualitatively suggests that the propellant may be expected to show 
about a five percent increase in modulus per year. Thus, for a maximum 
storage of 5 years the dynamic modulus may be expected to increase on 
the order of 25 percent. 

7.4 LIMITS OF PROPELLANT VARIABILITY 

The limits of expected propellant variability are established 
based on the analysis of the previous subsections. The master dynamic 
data for the 6 batches of propellant tested by CDS have been curve fit 
to the SRM propellant dynamic response model developed in [4] over 
the frequency range of 10 to 100 Hz with the results 
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G'(u>a T ) = 253 (ua T ) 0 ' 50 (r 2 = 0.999) (64) 

G"(wa T ) = 200 (^J 0,38 (r 2 = 0.97) (65) 

and 

tan6 = |r = 0.79 (ua T )"°* 12 (66) 

Figures 30 through 32 present the mean and + 3o standard deviation 
curves for G‘, G" and the loss tangent G"/G' , based on the analysis of 
variance results summarized in Table LI. 

Table HI summarizes the anticipated limits of the dynamic 
response of the SRM propellant at several frequencies. The + 3a 
limits estimated for tan6 are based on the + 3a curves for G' and G", 
since experience indicates that G" increases as G' does and vice versa. 
From Figures 30 and 31, and Table HI, it may be seen that G' may vary 
by as much as a factor of two and G" may vary by a factor of three. 

G' and G" may be expected to each increase approximately 25% with aging, 
but tamS should remain approximately unchanged with aging. 
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Temperature Reduced Frequency, 
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FIGURE 31. Master Dynamic Modulus, G", Versus Temperature Reduced Frequency, Mean and 3 Sigma Deviations. 
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DYNAMIC RESPONSE 


3UENCY 


G" 

tan6 

= G"/G' 


Mean 

+3 a 

-3a 

Mean 

+3o 

(psi) 

(psi) 




339 

550 

0.49 

0.67 

0.72 

419 

661 

0.48 

0.63 

0.68 

480 

724 

0.47 

0.61 

0.65 

56C 

832 

0.46 

0.58 

0.62 

624 

912 

0.45 

0.56 

0.60 

812 

1120 

0.44 

0.52 

0.55 

1005 

1350 

0.43 

0.48 

0.51 

1150 

1510 

0.42 

0.46 

0.48 
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VIII. ANALYSIS OF SRM DYNAMIC RESPONSE 

Dynamic response analyses of the Shuttle SRM have been conducted 
to investigate the sensitivity of the SRM frequency response to propel- 
lant variability. In addition, a heat conduction analysis has been 
conducted to determine if significant temperature gradients may exist 
in the SRM at launch time. 

8.1 GENERAL DISCUSSION OF VIBRATION ANALYSES OF SOLID ROCKET MOTORS 

General surveys of vibration effects on solid rocket motors have 
been presented by IBM [1,3], Hufferd and Fitzgerald [14], Fitzgerald 
and Hufferd [13], Achenbach [36], and Baltrukonis [37]. An extensive 
bibliography, including abstracts, is presented in [4]. 

Compared to other structures, solid propellant rocket motors 
possess rather unique structural features. Whereas the mass of the 
propellant is very large compared to the mass of the thin-walled case, 
the propellant is generally compliant and contributes little to the 
stiffness of the composite structure, affecting the gross dynamic 
behavior of the motor mainly through its mass. Furthermore, since 
propellants are viscoelastic materials, the propellant provides 
considerable damping to the system. Consequently, transient effects 
attenuate very quickly and steady-state oscillations require a high 
energy input since a large amount of energy is dissipated. 

In recent years finite element models have been developed for 
the steady state vibration of solid propellant rocket motors [38-42]. 
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The finite element models lead to definition of a stiffness matrix, a 
mass matrix and a damping coefficient matrix for the discretized system. 
Lumped-mass and consistent-mass methods are used for defining the mass 
distribution. The lumped-mass method concentrates the element masses in 
a manner that maintains the location of the center of mass of the 
structure, whereas, in the consistent-mass approach the mass is repre- 
sented as it is actually distributed in the structure. In practice, 
better agreement has been obtained using the lumped-mass approach. This 
approach requires a large numb*** of coordinate points for accurate 
analysis of systems with primarily distributed masses; however, the 
mass matrix for the entire structural assemblage is diagonal and positive 
definite, thereby reducing computation times. On the other hand, the 
consistenc-mass approach requires excessive computational effort to 
obtain a desired degree of accuracy. 

Damping matrices can be derived analogous to those used for the 
mass and stiffness matrices of appropriate internal damping characteris- 
tics. Propellant damping is accounted for using the dynamic complex 
representation of material properties. 

In dynamic problems, Hamilton's variational principle is frequently 
used to derive Lagrange's equations of motion for the discrete system. 

In the absence of darning the system of equations to be solved has 
the form 


[M] {r ( t ) } + [K]{r(t)> = (R(t)} (67) 

where the mass matrix [M] is composed of the element masses, {R ( t ) } 
represents the vector of nodal point forces at time t, and (r(t)) 
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represents the vector of nodal point displacements at time t. The 
natural frequencies of the system are obtained by taking the nodal point 
force vector R(t) to be the null-vector. 

For damped systems it is necessary to employ nodal damping to 
obtain solutions for the damped structure. Several methods are avail- 
able for introducing damping for a propellant grain. An effective 
damping coefficient, [C^] may be defined by 

[C eff ] = IC1 (68) 

where [K"] represents the imaginary part of the complex stiffness 
matrix. Some structural codes, such as NASTRAN, employ a structural 
damping coefficient, 5 , defined by 

C - = 7 tan 5 (69) 

when* G" and 9' are, respectively, the imaginary and real part of the 
complex shear modulus. 

Mora directly, steady state dynamic viscoelastic analysis may be 
obtained from the corresponding elastic counterpart, employing the 
viscoelastic-elastic correspondence principle, in which all elastic 
material constants are replaced by their corresponding viscoelastic 
counterpart expressed as frequency dependent complex numbers, i.e., 

E -*> E*(iu>) * E'(w) + 1E M (u) 

(j -► G*(iu) * G*(w) + iG"(w) 

K - K*(iu) * K' (u) + iK"(w) 

v .-*■ v*(1u) * ♦ iv"(u) 
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where K denotes the bulk nk )ulus and v represents Poisson's ratio. 
Lames relations still hold among the complex moduli in the frequency 
plane so that 


v* ( iu ) 



E*(ioj) 

6K*(iw) 


(70) 


An elastic material is represented in this notation by a complex modulus 
with a real p*rt which is constant with frequency and imaginary part 
which is zero for all frequencies. The common assumption of a constant, 
elastic bulk modulus implies that the dynamic Poisson's ratio will 
contain a nonzero imaginary part. Equation (67), for the natural 
frequencies of ibration, then takes the form 


or 


[-„ 2 M + K' + iK"] {r*} = 0 


-<u 2 M + K' 


K" 



! 2 

;*T M-K' 




(71) 


(72) 


The above replacement leads to a system of equations with twice 
as many unknowns as the correspoMing elastic problem; however, the 
solutions may be obtained in a manner analogous to the static case. 
Obvious simplifications are noted if the problem is handled throughout 
in complex arithmetic. 
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Two general techniques are employed for solving the equations of 
motion for vibration loading. The first consists of direct step-by-step 
integration of the simultaneous differential equations. The response 
history is divided into a sequence of time increments of equal length 
and the equations of motion are formulated on an incremental basis. 

The motion confuted during each time increment is added to the cond^t 1 ' s 
at the beginning of the increment to obtain the conditions at the end. 
Thus, the response is evaluated step-by-step through the desired time 
range, starting with any giver, initial condition. The acceleration is 
assumed to vary linearly during the time increment, which leads to 
expressions for the change in acceleration (and velocity) in terms of 
the initial conditions and the change of displacement. 

The incremental equations of motion are solved by standard static 
analysis methods; for example. Gauss elimination. 

In the second method of solving the equations of motion, the 
equtions are first decoupled by transforming to normal (mode shape) 
coordinates which are solved independently mode by mode, and the modal 
results superposed to obtain the total response. In this method, known 
as mode superposition, the displacement vector is expressed as a linear 
combination of mode shapes. Substitution of this expression into the 
equations of motion results in a single uncoupled equation for the n th 
mode of the system. The solutions cf each modal response equation are 
obtained by any convenient procedure. 
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8.2 PREVIOUS VIBRATION ANALYSES OF THE SHUTTLE SRM 

Several vibration analyses of the space shuttle SRM have been 
conducted by IBM [1,3] and Thiokol/Wasatch Division [6]. Additional 
one-eighth scale dynamic model tests have beer conducted by the NASA/ 
Langley facility [5]. Analyses of the one-eighth scale models has 
beer carried out by Grunman Aerospace Corporation [2]. The results 
indicate that reasonably good agreement is obtained from simplified 
analyses of dynamic response with more sophisticated NASTRAN analyses 
and experimental results. 

IBM [1] conducted simplified and NASTRAN analysis for, among 
other things, longitudinal through-the-thickness shear vibrations. They 
also investigated the effects of moduli and grain length. The properties 
;hown below, representative of the shuttle SRM, were used in their 
analysis for longitudinal shear vibrations: 

a = 20 inches 
b * 70 inches 
p p = 0.064 lb/in 3 

v p = 0.495 

They assumed the outer boundary to be rigidly clamped. In addition, 
based on information received from CSD, two different expressions were 
used for representation of the dynamic complex shear modulus, depending 
upon strain level. For low strains the representation' 

G 1 = 2800 uT 220 


was used in the’r simplified analysis, and 
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6* = 6244 (1 + 0.384i } 

was used In their NASTRAN complex eigenvalue analysis. A value of 
G' * 6244 was used In a NASTRAN real eigenvalue analysis. 

At hlqh strains, the corresponding representations were 

G' * 900 u>* 145 
and 

G* = 1365 (1 + .2451 ) 

The simplified analyses were carried out using the expression 

“» = 9 (n) 

where is a circular frequency coefficient which depends on the ratio 
a/b. Table LI 1 1 p«*esents the results for the first natural frequency 
using Equation (73) and a NASTRAN real and complex eigenvalue analysis, 
including the effects of length as determined from the NASTRAN analysis. 

IBM subsequently conducted a more complete NASTRAN analysis of 
longitudinal shear vibration including coupling with the case [3] and 
forced vibrations. These analyses treated several geometries which 
included more .realistic representations of propellant grain configura- 
tion with forward and aft skirts, nose cone and nozzle attached. A 
constant shear modulus of 1333 psi and a loss tangent of 0.246 were 
used in all analyses. The results indicated several longitudinal modes 
below 50 Hz with the first in the neighborhood of 15 Hz with a maximum 


shear strain of 0.95%. 
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TABLE LI II. RESULTS FOR THE FIRST NATIONAL FREQUENCY OF 

SHUTTLE SRH [3] 

Low Strains High Strains 


Equation (73) 

38.1 Hz 

17.9 Hz 

Real Modulus 

38.2 Hz 

17.9 Hz 

Complex Modulus 

38.8 Hz 

18.0 Hz 

L = 100 inches 

37.0 Hz 

17.3 Hz 

L * 150 inches 

37.7 Hz 

17.6 Hz 

L * 300 inches 

38.2 Hz 

17.9 Hz 


The NASA/Langley Research Center has conducted dynamic tests of 
one-eighth scale models of the shuttle SRM and mated external tank model. 
Three propellant grain configurations were manufactured by CSD using 
inert UTI-610 propellant which has the same binder-fuel -curative 
components as UTP-3001 propellant used in the TITAN III-C with inert 
sodium chloride and inert ammonium sulphate substituted for the live 
oxidizer, ammonium perchlorate. All had 0.1875 inch thick aluminum 
shells and were 19.5 inches in diameter and 147.32 inches long. The 
propellant length was approximately 145.4 inches. The propellant grain 
inner diameters were varied to represent lift-off. mid-burn and end-of- 
burn configurations. 

The model test at NASA/Langley [5] for the lift-off configuration 
indicated the first bending mode at 54.1 Hz, the first torsional mode 
at 135.3 Hz and the first longitudinal mode at 149.7 Hz. 
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Grumman [2] has also conducted dynamic analysis of the SRM 
models tested at the NASA/Langley facility usirg NASTRAN. They used 
a constant modulus of elasticity of E = 25,000 psi (G * 8333 psi) and 
a loss tangent of 0.52 based on material property data supplied by CSO. 

The results of Grumman 's NASTRAN analysis for the lift-off 
configuration gave the first bending mode at 56.2 Hz (compared to 58.4 
Hz for simple beam theory), 168.3 Hz for the first torsional mode 
(compared to 161 for single beam theory) and 196 Hz for the first 
longitudinal mode (compared to 180 for simple beam theory). 

Part of the reason for the discrepancy between Grumman' s analyses 
and the experimental results of NASA/Langley may be attributed to the 
fact that the material properties used by Grumman are probaDly unrealis- 
tic for the actual test conditions. 

The final set of analyses conducted on the shuttle SRM, to be 
discussed herein, were conducted by Thiokol /Wasatch Division [6] on a 
model 312 inches long, a grain 0,D, of 146 inches and an inner port 
diameter of 67 inches using properties of TP-H1123 and H-13 (inert 
TR-H1123) propellant, which are similar to the shuttle SRM prope'lant. 
Their results indicated a first longitudinal thickness shear mode at 
14.6 Hz with G' * 650 psi for TP-H123 propellant and 11.4 Hz with G' * 380 
psi for H-13 inert propellant at 7.0°F. The natural frequences varied 
from 8.6 Hz at 40°F to 26 Hz at 90°F for the first longitudinal shear 


mode. 
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8.3 SHUTTLE SRM HEAT CONDUCTION ANALYSIS 

Biot's approximation method [43] has been applied to cooldown 
of the Shuttle SRM in order to estimate the temperature distribution 
that would exist in the SRM at launch time. The analysis is axisymmetric 
so that the temperature. T. must satisfy the one-dimensional (radial) 
heat conduction equation: 


k 



pc 


3T 

at 


(74) 


where 



T 

3 

temperature 

r 

3 

direction of heat flux 

t 

3 

time 

c 

3 

specific heat capacity 

P 

= 

density 

k 

3 

coefficient of thermal conductivity 


Defining a heat flow vector H by 

k 3 2 T - nr II = 

7T ' p at ' arat 


it follows that 


or 



(75) 


H 


pcT dr 


(76) 



where q-j is the depth to which heat has penetrated at time, t. 
Employing variational techniques, one can define functions 


and 



(78) 


where V is the thermal potential, 0 is the dissipation, and Q is the 
thermal force which satisfy Lagranges equation 


av , 3D 

3qi 3q-j 


Q 


(79) 


Solution of Equation (79) yields the penetration depth, q^ , as a function 
of time, and hence the temperature spatial distribution as a function 
of time. 

The solution to Equation (79) is obtained for a circular port 
g**ain with inner radius a and outer radius b. The inner port is assumed 
insulated and a step temperature change is applied to the outer 
boundary. Two time regimes are considered: 
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(1) 0 < t < t. 

— a 

(2) t a <t<« 

where the time t. is the time for the heat front, defined as the distance 

d 

q.| from the outside, to reach the inner port. The input temperature is 
written as 


T(b,t) = Tj u(t) 


(80) 


where u(t) is the Heaviside unit step function. The spatial temperature 
distribution, based on experimental evidence and the form of the exact 
solution, is represented by 


T(r,t) 



for b-r<q 1 and 0<t<t a 
for b-r>q^ and t a <t<0 


(81) 


Substituting into Equations (75) through (78), Equation (79) becomes 




2b 


q l q l 


- 10 <b 


(82) 


for the penetration depth, q-j , where < = k/pc is the thermal diffus 
ivity. Equation (82) may be integrated directly to yield the cubic 
relation 


( 83 ) 
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Equation (83) has only one real root 

^1 b [ 2 cos (^ T ~) + ] 1 


(84) 


where 


cos 0 = 1 - 15 — x- 
b* 


Defining the thermal diffusion parameter, n 


2 


by 


2 

n 



Equation (84) may be written 


q 1 = b^ 11.03 n 2 = ™ nb = ^V<t 


(85) 


( 86 ) 


Setting q^ = b - a and defining X * b/a. Equation (86) may be solved 
for the time t a for the heat to reach the inner port, 

Q 



/X-l \2 
< ' X ' 


(87) 


or somewhat more accurately 




b 2 


2 A* 


- 3 X* 
“T 




In the limit as a -*■ 0, X ■+• ® for a solid cylinder and t, 

a 

For a thin cylinder, as a b, X -*> 1 and t. 0. 

a 



w • 


( 88 ) 
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In terms of the thermal diffusion parameter, n , Equation (88) 
may be written 


2 . 1 - 


2 

3X^ + 


ZX' 


3C X' 


(89) 


Once the temperature field reaches the inner port a new function 
must be defined in lieu of Equation (81) since the temperature at the 
inner port, T , must gradually increase. Since the inner port was 

a 

assumed insulated, the derivative of the temperature with respect to r 

2 2 

must be zero at r = a. Thus, for t > t,, or n > tv. the following 

— a — a 

temperature ^ield is appropriate 

T * T o(l-^K’-(& )2 ) + \ <*» 

Equation (90) may be solved in the same fashion as that used to 
determine q-j yielding 

r -0.214(f) 2 -Dl 

T a = T 0 h-e J (91) 


where 


-2 _ <t 

b n a 


2 

IL_ 


= t/t. 


(92) 


An inspection of Equation (91) reveals that it solves the appropriate 

2 2-2 

boundary conditions; that is, for n * n, , n =1 and T a * 0, and for 

a a 

n 2 - », T a - T 0 . 
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The nominal properties of the center segments of the shuttle SRM 

are: 

b = 72 in. 
a = 30 in. 

X = b/a = 2.4 
k 5 1 in^/hr 


Thus, assuming a step input temperature from cure at a nominal 140°F 
to 70°F at the case with the inner port insulated, the time for the 
temperature change to be felt at the inner bore is, from Equation (88), 
t = 6 days. From Equations (91) and (92), the ratio t/t, is given by 

a a 


t/t a = 1 - 4.67 

a 




Thus, the time for the inner port to reach the near equilibrium 
condition of 95 percent of the temperature change at the outside 
(i.e., T a = 74°F) , is 

O 


t 95 = 15 t a = 90 days * 

Thus, since the shuttle is not scheduled to be launched within 90 days 
of the manufacture of the SRM, it may be safely assumed that the SRM 
will be at a nearly equilibrium temperature equal to the local ambient 
temperature with only slight gradients through the propellant grain. 
These conditions represent the maximum cooling time, since internal and 
external cooling would achieve equilibrium more rapidly. 
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8.4 SRM DYNAMIC ANALYSIS 

Nine dynamic analyses were conducted following the procedures 
outlined in Section 8.1. Since temperature gradients of any consequence 
are not expected in the SRM, analyses were conducted at 60, 70, and 
80°F for the mean and + 3a variations in dynamic modulus. The 
pertinent analysis variables are shown in Table LIV. 

The analysis results and the corresponding propellant dynamic 
moduli are presented in Table LV. The results substantiate that the 
propellant contributes to the overall mass of the system, but does not 


TABLE LIV. DYNAMIC ANALYSIS VARIABLES 


a = 30 in. 
b = 72 in. 
t c = 0.5 in. 

L = 320 in. 

P c = 0.28 lb/in 3 

P p = 0.064 lb/in 3 

, E. = 29 x 10 6 psi 

is L_ 


TABLE LV. DYNAMIC ANALYSIS RESULTS, FIRST RESONANT 
FREQUENCY 


-3a 

Mean 

+3o 

Temperature 

(°F) 

warn 


03 

(Hz) 

(psi ) 

tan6 

03 

(Hz) 

a 

tanfi 

03 

(Hz) 

60° F 

1683 

0,424 

19.04 

2128 

0.48 

19.05 

2667 

0.50 

19.07 

70°F 

881 

0.455 

19-06 

1193 

0.56 

19 08 

1622 

0.591 

19.10 

80° F 

472 


19.10 

684 

0.628 

19.13 

1005 

0.672 

19.15 
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effect the overall stiffness. For moduli varying between 472 psi and 

pll 

2667 psi and structure damping coefficient, ; * — varying between 
0.21 and 0.34, the differences in the resonant frequency are indis- 
tinguishable. 

As a point of comparison, the simplified approximate analyses 
outlined in 4] were also conducted yielding resonant frequencies vary- 
ing between 10 Hz at 60°F and 26 Hz at 80°F. These values are in 
agreement with those obtained by Thiokol [6], but show more variation 
since a rigid casing is assumed in their determination. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 

Six batches of TP-H1148 propellant have been characterized for 
dynamic shear modulus and the influence of expected propellant variability 
on the dynamic response of the Space Shuttle SRM has been examined. 

Static and dynamic physical property characterization tests were 
conducted by the Chemical Systems Division of United Technologies 
Corporation on six batches of TP-H1148 propellant (TP-H1148-6,-7,-8, 
-9,9970096 and 9970115) furnished by Thiokol /Wasatch Division. The 
static tests consisted of constant elongation rate tests to failure and 
stress relaxation modulus tests at 40, 70 and 90°F. These tests were 
conducted to establish that the propellant batches tested were repre- 
sentative of propellant tested by Thiokol and to be used in the Shuttle 
SRM. Although mo**e variation was evident than in similar test results 
obtained by Thiokol, nominal results are the same jnd typical of PBAN 
propellants. 

Dynamic shear modulus tests were conducted at rious superimposed 
static (compressive) strain levels and several dynamic strain amplitudes. 
Typical of other solid propellants material nonlinearities were observed. 
The test data show stiffening of the dynamic modulus with increasing 
compressive static strain due to the reinforcing effect of the Ammonium 
Perchlorate and aluminum particles, and a softening with increasing 
dynamic shear strain level due to inherent norlinearities and debonding 
between the filler particles and the propellant binder. 
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The dynamic tests were conducted at static strain levels approxi- 
mating those expected in the SRM due to cure shrinkage, thermal cooldown 
grain slump, and ignition in pressurization. 

The test data generated by CSD and previously obtained data at 
Thiokol and the University of Utah on TP-HU48 and related PBmN 
propenants were analyzed in order to estimate expected batch-to-batch 
and within-batch propellant variability. Analyses of variance and a 
probabilistic analysis of the SRM propellant dynamic response indicate 
an expected one-sigma coefficient of variation of approximately 10 to 
15 percent. This variation yields a possible six- Id variation in 
dynamic modulus, however, from the dynamic response analysis of the SRM 
it is noted that the propellant provides so little stiffness to the 
overall structure that even a ten-fold charge in dynamic properties 
has no effect on the first resonance frequency. 

Since the space shuttle could be launched between 90 days and 
5 years following manufacture of the SRM, a heat conduction analysis of 
the SRM was conducted to determine the length of time required for the 
SRM to reach equilibrium ambient temperatures following cure. We also 
wanted to determine the effect a temperature gradient through the grain, 
if one existed at the time of launch, would have on propellant vari- 
ability and the SRM dynamic response. 3ased on the heat conduction 
analysis conducted it is estimated that the SRM will reach thermal 
equilibrium approximately 90 days after the completion of cure. There- 
fore, it is anticipated the temperature through the grain will be near 
equilibrium, perhaps oscillating slightly about the ambient temperature. 
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Based on the results of the thermal analysis, dynamic response 
analyses of the Shuttle SRM were conducted at (.0, 70 and 80°F using the 
expected + 3 and mean dynamic modulus values determined from the 
analysis of variance. These properties encompass a six-fold increase 
in dynamic modulus from low to high and a 50 percent change in the 
structural damping coefficient; but absolutely no differenc, in resonant 
frequency could be detected. Thus, in agreement with other studies, it 
is recommended that future dynamic analyses of the SRM consider only 
the stiffness of the case, but include the mass of the propellant. 

The results of the SRM dynamic analysis indicate that it is not 
necessary that any dynamic testing oe performed for batch acceptance. 

The normal static physical property and burn rate testing for quality 
control are sufficient to assure batch-to-batch reproducibility. It 
is thus recommended that dynamic testing not be conducted on a routine 
basis. 

The only dynamic loading environment not considered in this 
program or in previous programs is that arising from combustion 
instabilities. The most common instability in large solid rocket motors 
such as the Shuttle SRM is acoustic instability. The prevailing 
frequency modes are determined by the geometry and properties of the 
gas-and solid-filled rocket chamber. Acoustic mode vibrations in the 
gas and normal pressure fluctuations in the chamber cause a dynamic 
loading to the propellant. If the oscillations interact so as to cause 
resonance conditions in the propellant more pronounced pressure oscil- 
lations result, leading to an instability in burning and potential 
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over pressurization of the motor case. Since combustion instability 
is predominant in large motors under steady-state thrust conditions it 
is recommended that a study be undertaken to determine if combustion 
Instability represents a potential mode of structural failure for 
the SRM. 
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1.0 INTRODUCTION 

Under contractual agreement No. 76-001 CSD undertook the task of 
evaluating the dynaaic shear, constant rate to failure, and stress relaxation 
response for six batches of Space Shuttle SRM solid propellant (TP-H1 148) . 

These data are intended to be used to assess the effects of propellant vari- 
ability on the dynaaic response of the Space Shuttle SRM and to deteraine the 
accuracy with which the dynaaic response properties aay be known. 

The propellant was furnished to CSD in two shipeents totalling 234 lb. The 
first shipment was received 15 July 1977 on bill of lading No. B09063; the bal- 
ance was received on 17 November 1977 on bill of lading No. B11968. 
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2.0 WORK TO BE ACCOMPLISHED 

The speoifio work to be accomplished by CSD consisted of the following: 

A. Dynamic Shear Response 

Measure the dynamic shear modulus of six SRM solid propellant batches at 
various static and dynamic strain conditions. 

B. Constant Rate Failure Properties 

Measure constant rate properties of six SRM solid propellant batches at 2 
in. /min at 0°, 40°, 70°, 90°, and 120°F. 


C. Viscoelastic Relaxation Modulus 

Measure viscoelastic modulus of six SRM solid propellant batches at 40°, 
70°, and 90°F. 

All dynamic shear, constant rate, and stress relaxation tests have been 
completed and the test data are presented and analyzed in this final report. 
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3.0 TECHNICAL DATA 

The propellent evaluated during this oontraot consisted of material from 
one 600-gal hatch (TP-H1 148-9) and five ether batches. Their batch identifies' 
tions were -6, -7, -8, -9, -9970096, and -9970115. Procedures used for the 
evaluation are presented with data for each test mode. 

3.1 CONSTANT RATE DATA 

Constant displacement rate uniaxial tests were conducted using ICRPG class 
B samples ( ^ with plastic extensometers as shown in figure 3-1. The teats were 
conducted using an Instrou, with the standard strip chart recorder (figure 3-2) 
modified for automatic data reduction. The computerized output is presented 
as Eq, o b , e B , 0 C B , € \, and ^ where 

E * initial modulus 
o 

p 

o ■ r— (engineering stress) 
m A 

O 

e - strain at o 

> m 

a C » true stress (corrected, 7— (1 + £) 
m Am 

o 

c c 

e - strain at 0 

a m 

^ - strain at rupture 


The results for individual samples as well as mean and standard deviations 
of a group are listed on the computer printout. Additional details of the 
automated data acquisition system are presented in appendix A. 


ICRPG constant rate testa were conducted at 2 in. /min at 0°, 40°, 70°, 

90°, and 120°P on all six batches identified above. The failure properties are 
shown in figures 3-3 through 3-8. Five replloate samples tested per teat con- 
dition showed good reproducibility withir batches and reasonable batch-to-batch 
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Figure 3-2. Large Instron Tester with Hang Data Acquisition System 


variability. Nominal values of rupture strain, maximum stress, and initial 
modulus at T0°F are presen led in table 3-1 for all six batches. 


3.2 STRESS RELAXATION 

Stress relaxation tests were conducted at ^0°, 70° and 90°F using 1/2- by 
;1/2-in. bonded-end tensile specimens (figure 3-9) with a nominal 6-in. gage 
length. The test method described in section 2.2 of the 1CRPG Solid Propellant 
Mechanical Behavior Manual*' 2 was used for these tests. A strain of approxi- 
mately 3% was applied and stress decay with time was measured. The test strain 
levels. were determined from the gage length between the wood end- tabs; the 
increase in length was determined by optical cathetometer measurements (see 
figure 3-10). 
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Temperature. *F 

Figure 3-3. ICRPG Class B Data for Batch No. TP-H1 148-6 
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Temperature. *F 

Figure 3-4. ICRPG Class B Data for Batch No. TP-H 1148-7 
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Figure 3-6. ICRPG Class B Data for Batch No. TP-H1 148-9 
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Figure 3-7. ICRPG Claaa B Data for Bateh No. TP-H 1148-9970096 
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Temperature. *F 

Figure 3-8. ICRPG Class B Data for Batch No. TP-H1 148-99701 15 
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TABLE 3-1. UNIAXIAL CONSTANT RATE DATA FOR TP-H1148 PROPELLANT AT ?0°F 


T3177 



TP-H1148 Rupture Strain Maximum Stress (F/A ) Initial Modulus (Xa/z) 
Batch No. at ?0°F, Z at 70°F, psi at 70°F. psi 
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TABLE 3-2. STRESS RELAXATION TEST VALUES FOR INDIVIDUAL SAMPLES 

OF TP-H1I48 PROPELLANT 


Temperature 

40°F 70°F 90°F 


TP-H1148 


Strain, 

E ar 1 

Strain, 

E at 4 

Strain, 

E at 1 

Batch No. 


Z 

adn, psi 

Z 

min, psi 

Z 

min, psi 



2.88 

453 

3.03 

318 

2.97 

274 



? 03 

401 

2.47 

345 

2.76 

334 

6 

4 

3.04 

396 

2.58 

322 

2.99 

255 



2.85 

546 

2.59 

285 

2.80 

280 



3.C2 

512 

2.75 

285 

2.63 

274 

Average 


2.96 

462 

2.68 

311 

2.83 

283 



2.81 

634 

3.06 

374 

2.90 

444 



2.76 

606 

2.81 

422 

2.96 

475 

7 


2.69 

o42 

3.12 

370 

2.87 

377 



2.59 

4 6 

2.81 

392 

3.06 

332 



2.61 

J 

3.04 

382 

2.78 

456 

Average 


2.69 

615 

2.97 

388 

2.91 

417 



f 2.56 

462 

2.87 

375 

2.78 

282 



2.94 

400 

3.06 

313 

2.88 

265 

8 


I 2.8£ 

420 

2.98 

323 

2.80 

335 


i 

2.91 

493 

2.99 

278 

2.54 

340 



I 2.71 

504 

2.94 

312 

3.44 

228 

Average 


^ 2.80 

456 

2.97 

320 

2.89 

290 



r 








2.94 

563 

2.92 

404 

3.00 

323 



3.11 

404 

2.93 

338 

3.09 

299 

9 

i 

2.88 

409 

2.96 

419 

2.91 

323 



2.78 

624 

3.16 

283 

2.59 

297 



2.38 

524 

2.?6 

346 

3.06 

325 

Average 


2.82 

505 

2.99 

358 

2.93 

313 



3.09 

476 

3.04 

330 

3.16 

305 

9970096 

i 

2.89 

509 

3.02 

3^3 

2.74 

317 



. 2.62 

633 

3. 00 

314 

2.84 

412 



- 

- 

- 

- 

3.07 

304 

Average 


2.86 

539 

3.02 

326 

2.95 

334 



3.05 

563 

3.12 

379 

2.91 

362 

99701O 


2.98 

552 

3.19 

360 

2.91 

345 



. 2.95 

631 

3.02 

349 

2.94 

351 



- 

- 

- 

- 

2.78 

415 

Average 


2.5* 

58i 

3.11 

363 

2.88 

368 


Note: E • stress relaxation tdulus 
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Figure 3-11. Stress Relaxation Modulus at 3 % Strain for Batch No. TP-H1 148-6 
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Figure 3-12. Stress Relaxation Modulus at 3 % Strain for Batch No. TP-H1 148-7 
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Figure 3-13. Stress Relaxation Modulus at 3% Strain for Batch No. TP-H1 148-8 
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Figure 3-14. Stress Relaxation Modulus at 3* Strain for Batch No. TP-H1 148-9 
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Figure 3-16. Master Modulus Data for Batch No. TP-h i 148-99701 15 
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this shift did not adequately superimpose the data the shift factors were 
adjusted imperically . Both the ULF and iaperical shifts are shown where appli- 
cable. The aean curves froa the six batches are compared in figure 3-17 and 
the nominal 1-min moduli are given below: 


TP-H1148 Stress Relaxation Modulus 

Batch No. at 1 min t psi 


6 345 

7 460 

8 350 

9 395 

9970096 360 

9970115 385 


Sample-sample variance noted on the stress relaxation graphs reflects the 
within-carton gradient, potential sample machining and handling damage, aud 
experimental test uncertainties, A 6- in. gage length bonded-end sample was used 
to minimize the experimental errors caused by JANNAF dogbone flow through jaws 
and gage length errors . Actual strain measurements we~e obtained to 0.001 -in. 
accuracy using the optical cathetcmeter . Load calibre- cions were conducted 
using dead weight tests on each test channel. 

3.3 DYNAMIC SHEAR RESPONSE 

Small strain dynamic shear properties of solid propellants and other low 
modulus materials are routinely measured at CSD using a piezoelectric trans- 
ducer device. Illustrations of the dynamic test device and the associated 
electrical equipment are presented in figures 3«-l8 through 3-20. The device 
employs a stack of piezoelectric crystals as the dynamic driver source. This 
dynamic strain is transmitted through the solid propellant sample as shear and 
the output load is measured using a single piezoelectric crystal which func- 
tions as a load cell. 
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KM KM KK 10’ 10* 10» 10* KJ> 

Tun* (I/Ay), rmn 

Figure 3-17. Comparison of Stress Relaxation Modulus Curves for 


TP-H1148 Propellant 
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The complex aodulus of the test material can be represented for shear as 
either 


E* : E' + i B" 
or G» : G' ♦ i G" 


(3) 


In equation 3, E' is the real part and E" is the I mag in a r y part of the 
complex modulus. The complex modulus E* is given by the ratio of the maximum 
stress to the maximum strain. For low or medim damping materials E' is nearly 
equal to the elastic tensile modulus. The imaginary part of the modulus 
E" Is a damping term known as the loss modulus. The phase angle between 
the output force and driver crystal displacement is 5 and the loss tangent 
is defined by 

tan s or 2 ? (4) 

E* G' 
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For ocoposite propellants the visooelastlc modulus is generally 
represented by a Prony series 


E(t) 



(5) 


or G(t) • B q + e~ a i t (for shear) 

where the A's and B's are the Prony series coefficients and the e's are the 
time constants. 


The real and imaginary parts of the complex modulus for a linear visco- 
elastic material can be calculated from the Prony series using 

n A. 2 
E* (w) - Z 

i-0 w“+cx* 



However, this conversion does not work well for propellants which usually ex- 
hibit nonlinear modulus variations ^*5). 


The dynamics shear test device provides a direct measurement dynamic 
properties as a function of frequency and test temperature. The equipment is 
calibrated using stainless steel rings inserted between ^he piezoelectric 
driver stack (figure 3-19) and the readout crystal which are excited over a 
wide frequency range. During propellant tests the output signal and the phase 
angle are measured to determine the dynamic properties of the propellant over 
the desired frequency and test temperature ranges. 

The dynamic strain amplitude and the static compressive strain levels can 
be varied to determine the modulus sensitivity to various loading conditions. 

3.3*1 Dynamic Shear Propellant Data 

The dynamics shear modulus behavior was measured for each of the six 
batches of Shuttle SRM propellant. Production batch Ho. TP-H 1148-9 was 
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evaluated at three static and three dynamic strain levels. Test data at one 
strain level were obtained for each of the other batches of propellant at 90°, 
72° and 40°F. 

3.3.2 Batch Ho. TP-H1 148-9 

Dynamic modulus curves are presented in figures 3-21, 3-22, and 3-23 for 
bntdi No. TF Nil 48-9 at 0.001} dynamic strain and static strain levels of 3 }» 
6} and i2}. A - comparison of G' curves for each of these tests is presei ted in 
figure 3-24. Static strain sensitivity of batch No. TP-H1 148-9 at 100 Hz and 
72°F illustrates the dynamic modulus increase with compressive static strain 
levels: 


Compressive Strain, } G* , psl 

3 2,070 

5.8 2,950 

12 6,250 


where temperature = 72°P, and dynamic strain = 0.001}. 

Test data in figures 3-21, 3-22, and 3-23 shows that G* actually drops 
lower than G" at the 90°F test temperature for the 3} and 6} static strains. 
The soft propellant at 90°F is absorbing a larger percentage of the energy 
input to the sample or the sample may not be transmitting the load because of 
slippage. 

Empirical shift factors at the three temperatures are independent of 
strain levels (figure 3-25). The WLF equation is also the same since strain 
dependence is net considered in its derivation. 

Dynamic shear moduli as a function of dynamic strain level are presented 
in figures 3-26 through 3-29. The real part of the dynamic modulus did 
lnorease with decreasing strain level but the total change in modulus is less 
than 50}. This is a much smaller amplitude change than experienced with 
compressive static strain variatiors. 
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Figure 3-21. Master Dynamic Shear Modulus for Batch No. TP-H1 148-9 at 3} 

Static and 0.001 Dynamic Strains 
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Figure 3-22. Master Dynamic Shear Modulus for Batch No. TP-H1l4 r -) at 6} 

Static and 0.001} Dynamic Strains 
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Figure 3-23* Master Dynamic Shear Modulus for Batch Mo. TP-H1 148-9 at 12% 


Static and 0.001% Dynamic Strains 
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Figure 3-24. Comparison of Static Strain Le* sis for Batch No. TP-H1 148-9 
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Figure 3-26. 


Master Dynamic Sheer Modulus for Batch No. Tr-H1148-9 at 3% 


Static and 0.0001% Dynamic Strain 
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3-27. Comparison of Shear Modulus for Batch No. TP-H1148-9 at 2.9% 


Static and Different Dynamic Strain Levels 
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3-28. Comparison of Shear Modulus for Batch No. TF-H1148-9 at 5.8% 


Static and Different Dynamic Strain Levels 
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Figure 3-29. Comparison of Shear Modulus for Batch No. TP-H1148-9 at 12Z 
Static and Different Dynamic Strain Levels 
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Dynamio test results for the other batohes of propellant are pre sented in 
figures 3-30 through 3-3*1. These batohes are very similar to the major batoh 
(-9). The real part of the dynamio modulus for eaoh of the six batohes is 
presented in figure 3-35. 

The real part of the dynamio modulus at 100 Hz is given below: 

Batoh No. 0* , pal at 100 Hz 


9970115 

3*100 

-9 

2,120 

-7 

2,120 

-8 

2,100 

9970096 

1,900 

-6 

1,800 
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Figure 3-30. 


Frequency (cjAj). Hz 

Master Dynamic Shear Modulus for Batch No. TP-H1 148-6 
at 3X Static and 0.00 IX Dynamic Strains 
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Figure 3-31. 


Master Dynamic Shear Modulus for Batch No. TP-H1148-7 
at 3X Static and 0.00 IX Dynamic Strains 

14302 
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Figure 3-32. 


Frequency (u>Ay), Hz 

Master Dynamic Shear Modulus for Batch No. TP-H1148-8 
at 3% Static and 0.001% Dynamic Strains 
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Figure 3-33. Master Dynamic Shear Modulus for Batch No. TP-H1148-9970096 
at 3% Static and 0.001% Dynamic Strains 
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Figure 3-34. 


Frequency (cjAj). Hz 

Master Dynamic Shear Modulus for Batch No. TP-H1148-9970115 
at 3 X Static and 0.001Z Dynamic Strains 
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Figure 3-35. Betveen Batch Comparison of Propellant TP-H1148 
Real Part of Dynamic Shear Modulus at 3% Static and 
0.001Z Dynamic Strains 
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Batch No. 9970115 has the highest dynamio modulus with the others being close 
together; batch No. TP-H1 148-6 is the lowest. 

The real and imaginary part of the dynamic shear modulus values at 50 Hz 
are presented in table 3-3 for batch -9 and in table 3-4 for the other five 
batches. Individual values are listed along with mean values. 
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TABLE 3-3. REAL AND IMAGINARY DYNAMIC SHEAR MODULUS FOR BATCH NO. TP-H 1148-9 
• COMPARED AT 50 Hz 
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TABLE 3-4. SEAL AND IMAGINARY DYNAMIC SHEAR MODULUS FOR OTHER BATCHES OF TP-H1148 

PROPELLANT COMPARED AT 50 Hz 
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4.0 CONCLUSION 

Selected mechanical properties of six batches of 5RM solid propellant were 
measured using constant rate, stress relaxation, and dynamic shear test 
methods. Test data show very small batch-to-batch propellant variability and 
excellent sample-to-sample reproducibility. Stress relaxation modulus values 
are lower than the dynamic shear data, as expected for dynamic microstrain 
shear test conditions and the difference in time response. Dynamic propellant 
modulus values were sensitive to the static compressive strain level as well as 
dynamic strain. Some voids and fuel pockets were noted in the small dynamic 
shear test specimens which may explain the larger data variability than 
occurred with the larger test specimens used for the other test modes. 

The test strain levels covered 0.00 If to 3> and a correspondingly large 
modulus range. This strain range may be applicable to the Shuttle SRM dynamic 
loading conditions Specific modulus values for dynamic analysis can be 
selected from the available data by using the available NASTRAN computer 
results and defining the propellant strain range. 

The compressive strain condition of the dynamic shear sample resembles the 
SRM motor loading during ignition as the grain will be in compression from the 
interior ballistic and a shear load will be applied by gravity and accelera. 
tion. This particular dynamic shear test^®'^ is more appropriate than the 
Gottenburg disc and other dynamic tests where the strain fields experience 
large gradients through the sample and are not representative of SRh 
conditions. 
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1.0 AUTOMATED DATA ACQUISITION SYSTEM 

CSD currently Interprets Instron test data using an automated data acqui- 
sition system. This results in a low cost data processing system which p**o- 
vides on-line acquisition of test data which is more accurate and reliable 
than handreduced data. 

The equipment used in this automated data acquisition system is illus- 
trated in figure A-1. The Instron is mated with a CSD-developed electrical 
interface package which is designed to transfer load, strain, and time informa- 
tion into a Wang 700 series programmable calculator. The interface system 
conditions the analog and digital signals from the testing machine and formats 
the data for correct entry into the calculator. 

The test operator inserts the program into the Wang calculator using 
standard tape cassettes and manually inputs the test identification and test 
constants using the calculator keyboard. An output printer provides tabular 
and graphical test results and statistics in a report-type format. 

The Wang calculator-interface system is used to automate the ICRPG 
standard rate testing and data reduction. The test data (stress and strain) 
are collected, reduced, and presented as E 0 , <* m , * m , » c m , < c m , and « r . User 
options and calculation details are summarized below. 

1.1 DATA COLLECTION 

A Strain 

A constant crosshead rate (X) yields a simple relationship between strain 

(«) and elapsed time (T). 

« 3 XT /(gage length) (A-1) 


A- 2 



Tba system actually measures time on a 
60-cycle clock and calculates all strain 
using equation A-1 . 
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B . Stress 

The system oollects load information at 
uniformly spaced strain intervals (IC). 
The value collected is actually a per- 
centage (P) of a selected full-scale 
load range (R). The stress values are 
calculated by. 

« = (PR)/(samp.e area) (A-2) 

C. Strain Intervals 

Load information is collected at 
uniformly spaced strain intervals of 


Floor model Instron 



1.01* for X = 20 in./min U-3) 

6e- 0.60)1 for X = ? in./min 


Figure A-1. CSD-Automated 
JANNAF Data Acquisition 
System 


The calculator storage capacity is limited to a maximum of 13** load 
values which correspor to a maximum strain of 


e max s 135 * for X = 20 in./min (A-4) 

*max s 80J for x = 2 in./min 


These values seem reasonable for appropriate resolution over a wide 
range of test conditions. Unusual propellant properties or extreme 
test temperatures may find 6e too large (insufficient resolution) or 
too small (small t mM X or 6 load in the noise level). For X<5 in./min 
the operator can change 6e and consequently by inputting Y into 
register 26: 


y * nrr 


206 e 

X 



(A-5) 
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1.2 DATA REDUCTION 

A. Initial Modulus 

Three initial Modulus values are listed: E, El and E4. The calculation of 
these values involves a search through the stored stress data for the 
maximum Ac over a certain Ae range. The initial Modulus is defined 
as 


E 0 = 1.05 CAo/AeJaax (A-6) 

The three reported values differ in the A< (and corresponding Ao) range. E, 
El, and E4 are calculated over increments of Ae = 26, 6"» and 46, respec- 
tively. Since CSD's Method of discrete data collection se^aents the true 
stress strain curve, listing three values should give a better feel for 
the initial behavior of the curve. 

B. Zero Strain Point 

A zero strain point is necessary only to define other strains or at least 
one other strain to be used as reference. The zero strain point is 
indirectly located by assuming that the stress-strain curve, previous to 
the Initial modulus E, is a linear curve. From equation A-6 E and the 
stress value (a Q ) at the point where E is calculated are known, so the 
strain (e Q ) at this point is 

* 0 = o q /E (A-7 ) 

e Q is a reference point from which all other strains are generated. 

C. Maximum Stress and Strain 

The maximum stress value (o B ) is found by making a scan through tfc-t stored 
stress values. The strain at this point is found by counting the number 
of strain intervals (6«) from e 0 . Assuming M intervals, then, 

e m s e o * (A-8) 


A -4 
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Consequently, the ainiaua error associated with o B end i B Is 6 *. 

D. Maxiaua Corrected Stress end Streln 

The aaxiaua oorrected stress (o° n ) is defined es (O ♦ «)«1 ux ; it is 
found by soennlng for the largest product of stress x ( corresponding 
strain ♦ 1). «° m is defined es the strain at the aaxiaua corrected 
stress end found using equation A-8. The ainiaua error essooieted 
with o° m and «° n is s« . 

E. Failure Strain 

The failure strain (* P ) is defined as the strain at the rupture point, 
usually when o * 0. « r is found by using equation A-3 and is accurate 

to fit. 

1.3 OPTIONAL OUTPUT 

In addition to the standard output described above, two optional outputs 
are available. They are the following: 

A. Tabular Printout of «, E-j, o, and a° 

The segmented curve, consisting of the collected stress points between « Q 
and Intervals, with N<13*l (default) user inputed. For the first 
stress point to lie within eaeh Interval, the following is calculated and 
printed. 

e a strain at the first data pcint to lie ir. each interval (equation A-8) 

ET a (Ao /A«) with A« a 2 l< -and oentered above* 

o a value of first stress point to lie in each interval 
o° a ( 1 ♦ t ) ; using * and o froa above . 

The listed values will not neoessarily be equally increaented in 
strain (due to N-dependenoe) but rather the strain will correspond 
to actual stress data points (the first point in each interval). The 
last point is the failure point. 
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B. HD-20 Output 

The sepiented curve between < 0 and e p ia divided into 24 intervals and the 
flrat data point in each interval ia oonverted froa (e , a ) to (x, y) in., 
with ohart speed * X. The conversion ia 


X * «• (gage length) 

Y * lOo (area)/(load range) 


U-9) 


The 25th listed point is the failure point. The corresponding printout is 
in a format suitable for keypunching. 

The user operational keys for the Instron code are shown in table A-1 ; 
a typical printout ia shown in table A-2. 
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TABLE A-1 . USER OPERATIONAL KEYS OF IN ST RON CODE 
T3181 

Key Function Explanation Note 

01 Load A continuous loop is entered which Loop is exited 

calibration displays the load in the X-register by keying PRIME 

03 Print heading The heading is printed with appropriate 03 automatically 

program stops to allow manual typing calls 04 

04 Predata Gage length (2.7 in.}, XHD SPD (2 in./ 04 automatically 

min) and load range (20 lb) ars entered calls 08 
from the Wang keyboard 

08 Initialize Statistics are cleared for starting a 08 automatically 

series new series calls 10 

10 Test A program stop allows sample width 10 resets itself 

(0.375 in.) and thickness (0.5 in.) to for start of 

be entered; on keying GO, the test is a new test 

run and a summary printed 

12 Reject Rejects this sample from statistics 12 automatically 

calls 10 

13 Tabular A program stop allows entry of the max- 13 automatically 

printout imum number of data points to be printed calls 10 

out (default value of 134 is the maximum 
possible); this additional printout con- 
sists of e, E, S ando 

14 MD-20 Additional printout used as input to the automatically 

printout MD-20 code; the format is suitable for walls 10 

key-punching 

15 Statistics The mean and standard deviation for values 

of nonrejected samples is printed out; 
a program stop allows manual paper advance 
to a new page; key GO initiates summary 
printout starting with heading and leading 
into predata; keying 11 will skip the 
heading- predata printout and list only the 
samples 
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Th« 25 (x»y) pairs, described in the text, sre printed in keypunch forxat. 
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APPENDIX B 

DYNAMIC SHEAR CALCULATIONS AND COMPUTER OUTPUT 
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1.0 DYNAMIC SHEAR CALCULATIONS AND COMPUTER OUTPUT 

The dynamic shear modulus was computed for the real and imaginary parts of 
the oomplex modulus with the following equations: 

G * “ X & 008 ^ (Real) 

G" - B sin <j> (Imaginary) 

where L = sample thickness in inches (nominal 0.1 in.) 

A s sample area in square inches; contact area between propellant 

2 

and driver plate (nominal 0.18 in. ) 

p = calibration factor in lb/in. 

$ = phase angle between input and output signals. 

The P calibration curve versus output voltage is obtained from stainless 
steel rings cut from standard tubing with the following relationship 


F 




0,224 D 3 F/YLT 3 
F/AD - 4.46 YL(T/D) 3 
Young's modulus 
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The rings are made from standard tubing with values of 0.25 to 0.30 in., 
wall thickness of 0.001 to 0.065 in., and length of 0.050 in. This provides a vide 
range of load calibrations. 


The dynamic strain level is a linear function of the voltage applied to 
the driver stock of piezoelectric crystals and is calculated from 


2 KV 


(o-p) 1 


e = 


(B-3) 


where 


calibration constant 


L = sample thickness 


V (o p) = a PPl ied voltage 

e = dynamic strain output from driver crystal. 

A calibration curve is shown in figure B-1 . Dynamic shear data reduction 
for batch No. TP-1148-9 is shown in table B-1. 
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Figure B-1 . Calibration Curve 
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TABLE B-l. DYNAMIC SHEAR DATA REDUCTION FOR BATCH NO. IT- 1148-9 
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